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1.  Introduction 


In  Finsler  geometry,  each  point  on  the  base  manifold  can  be  envisioned  as  endowed 
with  a  vector  of  coordinates  denoting  its  position  from  the  origin  and  a  director  vec¬ 
tor,  also  referred  to  herein  as  an  internal  state  vector,  that  may  be  independent  of  po¬ 
sition.  Geometric  objects  such  as  metric  tensors,  connections,  and  derived  quantities 
(torsion,  curvature,  and  so  forth)  may  in  turn  depend  on  position  and  direction  or  in¬ 
ternal  state.  This  generality  is  in  contrast  to  classical  Riemannian  geometry,  wherein 
the  ultimate  dependence  of  such  objects  is  on  position  alone.  Indeed,  Finsler  geom¬ 
etry  encompasses  certain  geometries  of  Riemann,  Minkowski,  and  Weyl  as  special 
cases.  Its  generality  has  resulted  in  widespread  posited  field-theoretical  descriptions 
in  nearly  all  branches  of  physics:  general  relativity,1  gravitation,2  quantum  mechan¬ 
ics,3  electrodynamics,4  heat  conduction,5  and  the  mechanics  of  solids.6  The  latter 
topic  (i.e.,  continuum  mechanics/physics  of  deformable  bodies)  is  emphasized  in 
the  present  report.  In  this  context,  fields  describing  the  motion  of  material  particles 
comprising  a  body  must  be  introduced  along  with  evolution  of  the  internal  state, 
specifically  transformations  from  referential  or  Lagrangian  coordinates  and  initial 
state  vectors  to  spatial  or  Eulerian  coordinates  and  current  state  vectors. 

Finsler  geometry  is  attributed  by  name  to  the  doctoral  work  of  P  Finsler  nearly  a 
century  ago.7  Early  fundamental  contributions,  including  the  introduction  of  var¬ 
ious  connections,  were  set  forth  by  Cartan,8  Chern,9  and  Rund.10  Modern  mono¬ 
graphs  include  references  11-13.  Of  particular  interest  here  is  reference  14,  since 
it  includes  a  chapter  devoted  to  applications  in  finite  deformation  of  solids,  albeit 
with  content  limited  to  kinematics  alone.  See  also  the  historical  review  in  reference 
15 — a  paper  that  also  advances  Finsler  geometry  via  extension  of  the  Cartan-Clifton 
method  of  moving  frames — and  the  recent  categorization  of  Finsler  connections  in 
reference  16. 

Applications  of  Finsler  geometry  in  continuum  mechanics  and  physics  of  deformable 
solids  have  been  suggested,  but  not  fully  developed  or  realized,  since  the  middle  of 
the  20th  century.  Amari17  developed  what  appears  to  be  the  first  Finsler  geometric 
theory  of  deformation  of  solids,  applied  specifically  to  ferromagnetic  elastic-plastic 
crystals.  In  this  theory,  the  internal  state  vector  is  physically  linked  to  the  spin  di¬ 
rection  of  the  magnetic  moment,  and  dislocations  (a  fundamental  line  defect  in 
crystalline  solids1819)  are  associated  with  a  certain  torsion  tensor  related  to  anholo- 


Approved  for  public  release;  distribution  is  unlimited. 


1 


nomicity20  25  of  the  locally  relaxed  intermediate  state  of  the  crystal,  following  ear¬ 
lier  classical  differential-geometric  treatments  by  Kondo26  and  the  Japanese  school. 
Kondo27  briefly  discussed  possible  application  of  Finsler  geometry  to  describe  plas¬ 
tic  yielding.  Kroner28  and  Eringen29  suggested  how  Finsler  geometry  may  be  of 
potential  use  for  describing  mechanics  of  solids  in  the  context  of  generalized  con¬ 
tinuum  theories  such  as,  for  example,  reference  30,  but  did  not  further  develop  or 
expound  on  these  ideas.  Around  the  same  time,  Ikeda31,32  developed  a  theory  of  de¬ 
formable  media  with  close  connections  to  Finsler  space,  again  restricted  to  descrip¬ 
tion  of  kinematics  without  consideration  of  energy  functionals  or  equilibrium  equa¬ 
tions.  Apparently,  the  application  of  (pseudo-)Finsler  geometry  to  solid  mechanics 
remained  dormant  for  some  20  years  after  these  suggestions,  until  the  appearance  of 
work  by  Bejancu,14  followed  in  the  next  decade  by  contributions  from  Saczuk  and 
colleagues.6,33  (Herein,  a  space  is  designated  as  pseudo-Finslerian1416  rather  than 
strictly  Finslerian  when  a  fundamental  scalar  function  with  requisite  properties,11 
from  which  the  metric  tensor  is  obtained  by  differentiation,  does  not  exist.)  Theo¬ 
retical  developments  again  remained  scarce  for  15  years  following,  apart  from  some 
recent  work  on  anisotropic  acoustic  wave  propagation.34  For  a  more  comprehensive 
current  literature  review,  see  reference  35. 

It  is  speculated  that  Finsler  geometry,  in  contrast  to  Riemannian  geometry, 36-40  has 
heretofore  eluded  popularity  among  mechanicians  and  physicists  due  to  the  appar¬ 
ent  complexity  of  calculations,  despite  its  generality  and  descriptive  potential.  In¬ 
deed,  only  one  published  paper33  containing  solutions  to  a  boundary  value  problem 
in  Finsler-geometric  solid  mechanics  seems  to  exist,  and  these  solutions  were  ob¬ 
tained  numerically  rather  than  analytically  and  only  discussed  in  brief  (though  the 
somewhat  obscure  monograph41  contains  many  more  details  regarding  the  theory 
and  problems  considered  in  reference  33).  In  the  absence  of  solutions  to  physically 
meaningful  problems,  a  complex  new  theory  may  offer  little  advantage  or  insight 
over  simpler  existing  methods. 

The  purpose  of  this  report  is  to  initiate  a  new  theory  of  mechanics  of  deformable 
solids  with  microstructure  using  concepts  from  Finsler  geometry.  Although  a  few 
aspects  of  the  proposed  theory  are  drawn  from  prior  work,  notably  references  14 
and  33,  many  features  are  introduced  here  for  the  first  time.  The  theory  is  con¬ 
structed  with  an  aim  toward  obtaining  solutions  to  pertinent  boundary  value  prob¬ 
lems  in  mechanics,  physics,  and  materials  science.  Specifically,  problems  consid- 
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ered  herein  include  tensile  fracture  of  an  elastic  bar  (see,  for  example,  reference 
42  for  a  recent  analysis  in  a  different  multiscale  context),  slip  localization  in  an 
elastic  slab  under  simple  shear  (see,  for  example,  reference  43  for  a  recent  anal¬ 
ysis  via  phase  field  theory),  and  cavitation  (e.g.,  void  formation  and  expansion) 
in  a  spherical  elastic  domain.  The  fundamental  metric  tensor  entering  the  theory 
in  these  problems  accounts  for  microscopic  dilatation,  which  is  commonplace  in 
the  fracture  of  crystalline  rocks  and  minerals44-46  as  well  as  in  the  vicinity  of  dis¬ 
location  cores  in  crystals.18,47  To  reflect  such  local  volume  changes,  a  Weyl-type 
rescaling  (i.e.,  a  conformal  transformation)  of  the  metric  tensor  is  invoked.48,49  In  a 
novel  theory  of  thermal  stresses  based  on  Riemannian  geometry,  a  similar  rescaling 
of  a  metric  tensor  on  the  material  manifold  was  invoked  to  study  isotropic  non¬ 
linear  elastic  solids.50  Unlike  prior  theoretical  and  computational  studies  in  crystal 
inelasticity, 51-53  the  present  developments  do  not  require  a  multiplicative  decom¬ 
position  of  the  deformation  gradient  into  2  (or  perhaps  more)  terms,  but  such  a 
treatment  is  not  precluded  by  the  general  theory  and  has  been  proposed  elsewhere 
in  a  merging  of  continuum  phase  field  and  Finsler  geometric  treatments  of  defor¬ 
mation  twinning  in  crystals.35  Applications  considered  herein  are  also  restricted  to 
initially  homogeneous  bodies  (e.g.,  single  crystals  or  homogenized  polycrystals). 
Not  considered  explicitly  are  spatially  heterogeneous  bodies,  such  as  those  with 
varying  elastic  moduli51  or  varying  stress-free  strains,54  though  the  general  theory 
does  not  preclude  analysis  of  such  problems.  Though  this  report  does  not  apply  the 
theory  toward  problems  involving  phase  transitions  or  twinning,  as  often  studied  in 
the  context  of  martensite,55  the  theory  developed  herein  could  be  readily  applied  to 
such  problems  in  future  work. 

This  report  is  organized  as  follows.  Kinematics  and  geometry  pertinent  to  a  new 
pseudo-Finsler  theory  of  continuum  physics  are  described  in  Section  2,  including 
motions  at  macro-  and  microscales,  strain  metrics,  and  differential  geometric  ob¬ 
jects,  such  as  horizontal  and  vertical  connections.  An  energy  functional  over  the 
domain  is  developed  in  Section  3,  from  which  Euler-Lagrange  equations  yield  the 
conservation  laws  of  static  equilibrium  for  incremental  boundary  value  problems. 
Crucial  to  deriving  such  equations  is  an  extension  of  the  divergence  theorem  for¬ 
warded  in  reference  56.  Application  of  the  theory  toward  problems  involving  frac¬ 
ture  is  discussed  in  Section  4,  including  the  solution  of  a  1 -dimensional  (1-D)  tensile 
decohesion  problem.  Application  toward  problems  involving  slip  (e.g.,  dislocation 
glide,  adiabatic  shear  band  formation,  or  mode  II  separation)  is  presented  in  Section 
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5.  Application  toward  problems  involving  cavitation  is  given  in  Section  6.  In  the  ap¬ 
plications  sections,  links  among  geometric  objects/parameters  entering  the  theory 
and  physical  concepts  are  highlighted,  and  comparisons  with  predictions  of  other 
classical  field  approaches  are  given.  Conclusions  follow  in  Section  7. 

2.  Pseudo-Finsler  Geometry  and  Kinematics 

Discussed  in  turn  are  aspects  of  the  referential  (e.g.,  initial)  configuration  of  the 
material  body,  the  deformed  configuration  of  the  material  body,  and  then  the  trans¬ 
formations  (e.g.,  motions)  between  the  2  configurations. 

2.1  Reference  Configuration  Geometry 

The  reference  configuration  is  identified  with  a  particular  instant  in  time  at  which 
a  deformable  solid  body  is  considered  undeformed,  following  the  usual  conven¬ 
tion  of  continuum  physics.19  A  differential  manifold  971  of  spatial  dimension  3  is 
then  physically  identified  with  a  deformable  solid  body  embedded  in  ambient  Eu¬ 
clidean  3-space.  Let  A"  G  971  denote  a  material  point  or  material  particle,  and  let 
{XA}(A  =  1,  2,  3)  denote  a  coordinate  chart  that  could  be  assumed,  in  the  inter¬ 
est  of  brevity,  to  completely  cover  971.  The  body  manifold  may  be  taken  as  simply 
connected  herein  for  simplicity,  though  such  an  assumption  is  inessential;  in  fact, 
the  theory  is  later  applied  to  describe  breakage,  slip,  and  cavitation  processes  that 
may  preclude  simple  connectivity  of  the  body  in  one  or  more  conditional  states.  At¬ 
tached  to  each  point  is  a  vector  D,  or  equivalently,  a  chart  of  secondary  coordinates 
{Da}(A  =  1,  2,  3)  is  assigned  that  is  treated  as  a  field  description  of  microstruc¬ 
ture  in  the  solid  and  can  be  associated  with  a  second  manifold  if  of  dimension  3. 
Herein,  D  need  not  be  of  unit  length.  Analogoulsy,  one  may  assume  for  temporary 
convenience  that  if  is  simply  connected  and  covered  by  a  single  chart,  though  again 
this  assumption  is  not  essential.  Regarding  notation,  dependence  of  a  function  on 
(A",  D)  implies  dependence  on  charts  ({A^4},  { DA }). 

Following  the  notation  of  reference  14,  the  description  of  the  reference  state  of  the 
body  can  be  couched  in  terms  of  pseudo-Finsler  geometry.  Define  Z  =  (3,  II,  Oft,  if) 
as  a  fiber  bundle  of  total  (pseudo-Finsler)  space  3  (dimension  6),  where  II  :  3  — >■  971 
is  the  projection  and  it  the  fiber.  A  chart  covering  3  is  then  {A",  D}.  The  natural  or 
holonomic  basis  on  3  is  the  field  of  frames  {  ■  jjjja  }  •  Coordinate  transformations 
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from  {X,  D}  to  another  chart  {A",  I) }  are  of  the  Finsler  form11 


xA  =  xA(x\x2,  x3),  da{x)  =  q^(x)db.  (i) 


Let  Qr  =  jjYTj-  From  the  chain  rule,  holonomic  basis  vectors  on  Tj  then  transform 
as11-14 


d 


dXB  d 


d2XB 


D 


c 


d 


d  dXB  d 


dXA  3XAdXB  dXAdXc  d  DB'  dDA  dXAdDB 


(2) 


Let  N§(X,  D )  denote  nonlinear  connection  coefficients.  Nonholonomic  basis  vec¬ 
tors,  which  unlike  their  holonomic  counterparts  transform  as  typical  vectors,  are 


5  d  atB  d 
SXA  8XA  A  dDB  ’ 


5Da  =  dDA  +  NgdXB . 


(3) 


Noting  the  scalar  products  w,dXA )  =  5^  and  (t^b,SDa)  =  5^,  the  set 

{^Ya,  serves  as  a  convenient  local  basis  on  T3,  and  likewise  the  reciprocal 
set  {dXA}  SDa}  for  T*j.16  The  Sasaki  metric  tensor  invokes  the  latter: 


G(X,  D)  =  Gab(X,  D)dXA  ®  dXB  +  Gab{ X,  D)5Da  <g>  6DB.  (4) 


Components  GAB  and  their  inverse  components  GAB  are  used  to  lower  and  raise 
indices  in  the  usual  manner,  and  G(X,D)  =  detfG'  udA".  /))].  Differentiation  is 
hereafter  denoted  by  the  following  condensed  notation: 


9a(-) 


m.  a  ( )  =  M. 

dXA  ’  AK)  8DA ’ 


M‘) 


M. 

5XA 


dA(-)-NBdB(-). 


(5) 


Christoffel  symbols  of  the  second  kind  for  the  Levi-Civita  connection  on  5DT  are 
derived  in  the  usual  way: 


7 bc  —  \GAD (9cGbb  +  8bGcd  ~  dBGBc )  —  GADryBCD-  (6) 


Cartan’s  tensor  in  the  reference  configuration  is  defined  as 

CBC  =  \GAD{dcGBD  +  BbGcd  —  9dGBc )  =  GadCBcd-  (7) 
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Horizontal  coefficients  of  the  Chern-Rund  and  Cartan  connections  are  defined  as 


Tfic  —  \Gad{5cGbd  +  $bGcd  —  $dGbc)  —  GADV  bcd-  (8) 

Components  of  the  spray  and  derived  nonlinear  connection  coefficents  are,  respec¬ 
tively, 


=  \^cDBDc,  Ga  =  BbG/ 


(9) 


Letting  V  denote  the  covariant  derivative,  horizontal  gradients  of  basis  vectors  are 
determined  by  the  generic  affine  connection  coefficients  HBC  and  KBC : 


*  =Hi  6 


V7  d  KA  B 

Vs/6XbWT^  ~  KBC~: 


(10) 


s/sxB sxc  -  ±±bc~8Xa'  ' d 7 dX~  dDc  ~  ^BCd DA' 

Vertical  gradients  are  denoted  by  the  generic  connection  coefficients  VBC  and  YBC: 

(ID 


v  d  -  vA  9 

Vd/dDBgJ^  -  VBC 


dDA ’ 


SXA' 


Developments  to  this  point  apply  for  pseudo-Finsler  space  or  Finsler  space.  The 
latter  classification  holds  when  a  C°°  fundamental  scalar  function  £(X,  D )  exists 
at  every  point  of  it  \  0,  homogeneous  of  degree  one  in  D.1 1  from  which  the  metric 
tensor,  spray  connection  coefficients,  and  Cartan  tensor  are  derived,  the  latter  with 
additional  symmetry  not  necessarily  present  in  Eq.  7: 

Gab  =  \BaBb{£2)-,  Cabc  —  !<9a<9b(9c(£2); 

GAb  =  7 £CDC  ~  Cac1cdeDdDe  =  YabcDc. 

Formally,  the  Chern-Rund  connection  is  defined  when  Eq.  12  holds  and  NB  =  GB, 
Hbc  =  Kbc  =  TAC,  and  VBC  =  YBC  =  0;  the  Cartan  connection  is  defined  when 
Eq.  12  holds  and  NA  =  GB,  HAC  =  KAC  =  VAC,  and  VAC  =  YAC  =  C£c.16  Let 
(•)|c  denote  horizontal  covariant  differentiation  in  a  coordinate  chart  {A"6'}.  Then 
when  either  of  these  connections  is  used,  the  horizontal  covariant  derivative  of  the 
metric  tensor  vanishes: 

^ AB\C  ~  OC^AB  ~  1  CA^DB  ~  ±  CB^DA 

=  BcGab  —  Nb8dGab  —  r  caGdb  —  FbbGda  =  0. 

A  (pseudo)-Finsler  space  degenerates  to  a  Riemannian  space  when  Gab{ X,  D)  — y 
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Gab(X),  and  to  a  locally  Minkowskian  space  when  £(X,  D )  — >  £(D)}6 

Let  dX  denote  a  differential  line  element  on  Oft,  and  let  d D  denote  a  corresponding 
element  on  It,  both  with  components  referred  to  the  nonholonomic  basis.  Squared 
lengths  of  these  elements  with  respect  to  Eq.  4  are 

|dA|2  =  (dX,  GdX)  =  Gab dXAdXB ,  |d£>|2  =  (d D,  GdD)  =  GABdDAdDB. 

(14) 

The  traditional  scalar  volume  element  and  the  corresponding  volume  form  of  DJt 
are56 

dV  =  VGdX'dXMX3,  (in  =  VGdX1  A  dX2  A  dX 3.  (15) 

The  area  form  corresponding  to  a  compact  region  of  Wl  is 

Q  =  y/BdU1  A  dU2.  (16) 

The  embedding  of  dDJl  in  911  is  represented  by  the  local  parametric  equations  XA  = 
XA(Ua)  (a  =  1,2),  Ba  =  and  B  =  det {BaGabB%).  See  references  57 
and  58  for  a  comprehensive  treatment  of  surfaces  in  (Riemannian)  geometry.  The 
following  identities  are  also  noted: 

Min  VG)  =  tbABi  (Vg)\A  =  dA(VG)  -  nb8b(Vg)  -  Vghbb.  (17) 

Stokes’  theorem  in  terms  of  a  generic  C 1  differentiable  form  a  is 


Let  a(X,  D )  =  VA(X,  D)Q(X,  D )  be  a  2-form  in  Eq.  18,  and  let  VA  be  contravari- 
ant  components  of  vector  field  V  =  MMr-  Let  the  horizontal  connection  be  such 
that  the  second  of  Eq.  17  vanishes;  for  example,  HBB  =  TAB  (y/G)\A  =  0.  Then 
in  a  coordinate  chart  (A""4},  Eq.  18  becomes56 

[  [V\i  +  {vAccBC  +  8BvA)DBA}dn=  I  vAxAn,  (19) 

Jm  J  om 

where  N A  is  the  unit  outward  normal  to  the  domain  of  integration,  V^A  =  5AVA  + 
VaHba,  and  Dba  =  8ADB  +  NB. 
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2.2  Deformed  Configuration  Geometry 


The  current  configuration  is  identified  with  a  particular  instant  in  time  at  which  a 
solid  body  is  considered  deformed.  The  present  discussion  fully  parallels  that  of 
Section  2.1  but  with  an  adjustment  in  notation  for  deformed  coordinates  and  their 
indices,  which  are  denoted  via  lowercase  rather  than  capital  fonts.  A  differential 
manifold  m  of  spatial  dimension  3  is  identified  with  a  deformed  solid  body  em¬ 
bedded  in  ambient  Euclidean  3-space.  Let  x  G  m  denote  a  spatial  point,  and  let 
{xa}(a  =  1, 2,  3)  denote  a  coordinate  chart  that  is  assumed  to  completely  cover  m, 
which  may  correspondingly  be  taken,  by  inessential  assumption,  for  now  as  simply 
connected.  Any  assumption  of  simple  connectivity  will  be  relaxed  later  in  the  con¬ 
text  of  applications  involving  discontinuous  motions  associated  with  fracture,  shear 
banding,  and  cavity  formation;  in  such  cases  an  atlas  consisting  of  multiple  charts 
may  be  needed  to  adequately  cover  the  deformed  body.  Attached  to  each  point  is  a 
vector  d ,  or  equivalently,  a  chart  of  secondary  coordinates  {da}(a  =  1,  2,  3)  is  as¬ 
signed  that  is  treated  as  a  field  description  of  micro  structure  and  can  be  associated 
with  a  second  manifold  u  of  dimension  3.  Herein,  d  need  not  be  of  unit  length.  One 
may  now  likewise  assume  u  is  simply  connected  and  covered  by  a  single  chart  for 
simplicity,  though  the  theory  developed  in  this  report  may  be  readily  applied  toward 
bodies  that  are  not  simply  connected. 


Again  following  the  notation  of  reference  14,  the  deformed  state  of  the  body  can 
be  couched  in  terms  of  pseudo-Finsler  geometry.  Define  f  =  (3,  it,  m,  u)  as  a  fiber 
bundle  of  total  (pseudo-Finsler)  space  3  (dimension  6),  where  it  :  3  — >■  m  is  the 
projection  and  u  the  fiber.  A  chart  covering  3  is  {. x ,  d}.  The  natural/holonomic  basis 
on  3  is  {gf^,  Coordinate  transformations  from  {a;,  d)  to  another  chart  {x,  d} 
are 

xa  =  xa(x1,x2,xi),  da(x)  =  q%(X)db.  (20) 

Letting  qb  =  holonomic  basis  vectors  on  T3  transform  as 

d  dxb  d  d2xb  ~  d  d  dxb  d 

_ = _ | - dc _  _ _  = _ .  (21) 

dxa  dxa  dxb  dxadxc  ddb  ’  Qda  dxa  ddb 


Let  nb(x,  d)  denote  nonlinear  connection  coefficients.  Nonholonomic  basis  vectors 
are 


S  d  b  d 
Sxa  dxa  ,?a<9  db' 


5da  =  dda  +  nabdxb. 


(22) 


Approved  for  public  release;  distribution  is  unlimited. 


8 


The  set  {^,  ^ }  serves  as  a  convenient  local  basis  on  T3,  and  {dxa,  5da}  for  T*%. 
The  Sasaki  metric  tensor  in  spatial  coordinates  is 


g(x,  d)  =  gab(x,  d)dxa  <8>  dxb  +  gab(x,  d)5da  <8>  8db. 


(23) 


Components  gab  and  inverse  components  gah  are  used  to  lower  and  raise  indices, 
and  g{x,  d)  =  det [gab(x,  d)].  Spatial  differentiation  is  written  in  condensed  form  as 


Q  (.)  =  Q  (.)  = 

oU  dxa ,  a[)  dda’ 


W  =  S£  =  BM  -  »&(•)• 


(24) 


Christoffel  symbols  of  the  second  kind  for  the  Levi-Civita  connection  on  m  are 


7 be  =  yad(9cgbd  +  dbgcd  -  ddgbc)  =  gad^bcd- 


(25) 


Cartan’s  tensor  in  the  current  configuration  is  defined  as 


Cbc  =  w  \dcgbd  +  dbgcd  ~  ddgbc)  =  g  Cbcd. 


(26) 


Horizontal  coefficients  of  the  spatial  Chem-Rund  and  Cartan  connections  are 


r6C  =  v>c9bd  +  hged  -  Sdgbc)  =  g  r bcd. 


(27) 


Spatial  components  of  the  spray  and  derived  nonlinear  connection  coefficents  are, 
respectively, 


g°  =  Hcdbdc ,  9b  =  dbga.  (28) 

Letting  V  denote  the  covariant  derivative,  horizontal  gradients  of  basis  vectors  are 
determined  by  the  generic  affine  connection  coefficients  //,'(.  and  K'f'jC: 


^  d  TTCL  d  _  d  T.a  d 

^ 5 ^5xb~6xc  ~  ^bc^^a'  ^S/Sxb  =  d^bc 


'Sx° 


ddc  dda' 

Vertical  gradients  are  denoted  by  the  generic  connection  coefficients  Vbc  and  Ybac: 


(29) 


V  —  =  Va  — 

Vd/ddddc  Vbcdda' 


V  b  6  -Ya  6 


(30) 


The  spatial  configuration  correlates  with  Finsler  rather  than  pseudo-Finsler  space 
when  a  C°°  fundamental  scalar  function  f(x,  d)  exists  at  every  point  of  u  \  0,  homo- 
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geneous  of  degree  one  in  d,  from  which 


gab  =  ldadb(i2), 


gab  =  rbcdc  ~  Cabclcdeddde  =  r abcdc 


Cabc  =  \dadbdc{i2). 

(31) 


The  Chem-Rund  connection  is  invoked  when  Eq.  31  holds  and  nb  =  gb,Hbc  = 
Kbc  =  T£c,  Vbc  =  Ybc  =  0;  the  Cartan  connection  is  invoked  when  Eq.  31  holds  and 
nb  =  gb,  Hbc  =  Kbc  =  Ybc,  Vbc  =  Ybac  =  Cbc.  Let  (-)|c  denote  horizontal  covariant 
differentiation  in  a  spatial  chart  {a;c}.  Then  when  either  of  these  connections  is 
used,  the  horizontal  covariant  derivative  of  g  vanishes: 


gab\c  dcgab  Y  cagdb  Ycbgda  dcgab  iic  O^gaii  Ycagdb  Ycbgda  0.  (32) 


A  (pseudo)-Finsler  space  degenerates  to  a  Riemannian  space  when  gab(x,  d) 
gab(x),  and  to  a  locally  Minkowskian  space  when  l(x,  d)  — >  1(d). 

Let  dx  denote  a  differential  line  element  on  m  and  dd  denote  a  corresponding  ele¬ 
ment  on  u.  Squared  lengths  of  these  elements  with  respect  to  Eq.  23  are 

|cEc  | 2  =  (dx,gdjt)  =  gabdxadxb ,  |ck/|2  =  (dd.gdd)  =  gabddaddb.  (33) 


The  scalar  volume  element  and  volume  form  of  m  are 


dv  =  v/^da:1da;2da:3,  du  =  yfgdx1  A  da:2  A  da:3.  (34) 

By  a  simple  change  of  notation  from  referential  to  spatial  quantities,  an  area  form 
c o  can  be  introduced  analogously  to  Yt  in  Eq.  16,  as  can  spatial  versions  of  the 
coordinate-free  Stokes’  theorem  in  Eq.  18  and  Rund’s  horizontal  divergence  theo¬ 
rem  in  Eq.  19. 

2.3  Deformation  Kinematics 

Transformations  from  referential  to  spatial  coordinates  (59T  to  m)  and  vice  versa  are 
denoted  by  the  C 2  functions 

xa(X,  D)  =  ifa[X ,  D(X)\,  XA(x,  d)  =  d(x)(.  (35) 

Since  the  present  theory  is  quasi-static,  time  does  not  enter  such  functions  as  an 
explicit  independent  variable.  Incorporation  of  the  internal  state  ( D  or  d)  in  these 
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motion  functions  distinguishes  Finsler  kinematics14-33  from  the  usual  kinematics  in 
Riemannian  geometry  of  classical  continuum  physics.25,59  State  vector  mappings 
are  of  the  affine  form 


da(X,  D )  =  i)a [X,  D{X)}  =  daB (X)Db,  Da(x ,  d)  =  0A[x,  d(x)]  =  0A(x)db. 

(36) 

The  deformation  gradient  from  reference  to  current  (pseudo)-Finsler  tangent  spaces 
is  defined  as  the  delta  derivative 


F(X,  D)  =  FHX,  D)X  ®  ^  0  (ixA 


F%  =  SA<pa  =  dAxa  -  N*dBx 


Sx{X,  D ) 

8X  ’  (37) 


The  analogous  mapping  from  spatial  to  referential  tangent  spaces  is  the  2-point 
tensor 


f(x,  d)  =  fA(x,  d)J^ T  <g>  dxa  =  —  ^  <g>  dxa  =  SX^X’  ^ 


SXA  ~  5xa 

fA  =  8a$A  =  daXA  -  nbaBbXA. 


8XA 


5x 


(38) 


Herein,  for  a  2-point  tensor  such  as  F  in  Eq.  37  and  using  Eq.  35,  define 

»  pml  _  9x‘(X,D)  dx°(X,D)dDB(X) 
dAx  [X,D(X) ]  -  QXA  +  —555 - 

with  an  analogous  definition  of  daXA[x ,  d(x)]  in  Eq.  38: 


daXA[x,d(x)  = 


dXA(x,  d)  dXA(x,  d)  d db(x) 


dxa 


+ 


ddb 


dx° 


(39) 


(40) 


Transformation  equations  for  differential  line  elements  follow  by  generalizing  fun¬ 
damental  postulates  of  continuum  mechanics  in  Riemannian  space25-59  to  account 
for  the  nonholonomic  bases  of  Finsler  space,  whereby  partial  coordinate  derivatives 
are  replaced  with  delta-derivatives: 

r  c 

dx  =  ^dX  dxa  =  F% dXA,  dX  =  ^dx  ^  dXA  =  fAdxa.  (41) 

It  follows  that  volume  elements  and  volume  forms  transform  between  reference  and 
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spatial  frames  as 


du  =  JdV  =  [det{FX)y/g/G\dV,  doo  =  JcM; 
dV  =  jdv  =  [det(fA)  y/i G/g]dv ,  dQ  =  jduj. 

Lengths  of  deformed  and  initial  line  elements  can  be  compared  using  the  deforma¬ 
tion  tensor  C: 

| cLr  | 2  =  F%FbBgab dXAdXB  =  CAb dXAdXB  =  (dX,  CdX),  (43) 

where 

C  =  CABdXA  ®  dXB  =  F1gabFbBdXA  ®  dXB .  (44) 

It  follows  that  clet (6(4 R)  =  J2G.  For  the  directors  or  state  vectors,  a  similar  con¬ 
struction  using  Eq.  36  gives 

\d\2  =  rda  =  E  abDaDb  =  (. D,ED ),  (45) 


with 

S  =  Eab5Da  ®  5Db  =  rAgabdbB5DA  ®  SDB.  (46) 

A  transformation  rule  for  gradients  of  nonholonomic  bases  is  obtained  from  defini¬ 
tions  Eq.  29  and  Eq.  37: 


S  5xa  S  S  h  S 

Vs/SXAS&  =  6X AV5/SXa8^  =  FAW&/Sxafa~c  =  FAHacJ^b- 


(47) 


This  transformation,  which  is  used  later  in  the  total  covariant  derivative  operation25 
extended  here  to  2-point  tensors  in  Finsler  space,  provides  further  motivation  for  the 
presently  proposed  definition  of  the  deformation  gradient  F  as  a  delta-derivative. 


3.  Energy  Functional  and  Conservation  Laws 

The  following  variational  principle  is  set  forth,  where  T  is  the  action  integral  for  a 
closed  and  simply  connected  region  of  9JT  with  boundary  f/JJl,  and  surface  forces 
are  p  =  padxa,  a  mechanical  load  vector  (force  per  unit  reference  area),  and  z  = 
zaSDa,  a  thermodynamic  force  conjugate  to  the  internal  state  vector: 

6*(<p,D)=<[  {(p,6x)  +  (z,SD))n.  (48) 

JdW 
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Letting  i/j  denote  potential  energy  density  (i.e.,  internal  or  free  energy  density  in  the 
absence  of  kinetic  and  thermal  effects),  this  becomes 

5  f  =  < J  [pa8xa  +  zA8(DA)]Q,  (49) 

Jm  Jam 

where  the  first  variation  of  D  is  enclosed  in  parentheses  to  avoid  confusion  with 
basis  vector  8DA.  The  following  functional  form  of  the  energy  density  per  unit 
reference  volume  on  9JT  is  assumed: 


VD,  G)  =  Da ,  Da,  Gab).  (50) 


This  form  of  the  energy  density  is  motivated  by  phase  field  theory,43-60  whereby 
the  internal  state  vector  D  is  treated  here  analogously  to  an  order  parameter.  Fur¬ 
ther  motivation  is  obtained  from  continuum  mechanical  models  of  liquid  crystals, 
wherein  director  vector  gradients  may  enter  the  thermodynamic  potentials.61  Later 
in  Sections  4,  5,  and  6,  a  more  precise  physical  meaning  will  be  assigned  to  the  in¬ 
ternal  state  vector  in  the  context  of  example  problems.  Thermodynamic  forces  are 
introduced  by  taking  the  first  variation  of  Eq.  50: 


*+ = + w*S(dA)  +  + ^LSGab 


(51) 


=  PA8FaA  +  Qa8(Da )  +  Zba8D^b  +  S  8GAb- 

Imposition  of  spatial  coordinate  invariance  leads  to  the  restricted  form 


=  i!>[C{F,g),D,VD,C\  =  iP(Cab,Da,Da,Gab),  (52) 


from  which  the  first  Piola-Kirchhoff  stress  PA  and  Cauchy  stress  aab  obey 

P?  =  29o i.FbB^P,  =  2jF‘AFl PF  =  (53) 

Symmetry  of  Cauchy  stress  is  consistent  with  the  balance  of  angular  momentum  of 
classical  continuum  mechanics.19'59 


Noting  that  variation  <5(-)  is  performed  with  X  fixed  but  D  variable, 


5F%  =  8A(8ipa)  -  dB<padcN*8(Dc), 


(54) 
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(55) 


SDfB  =  [S(Da)]]b  -  ( BCNA  -  BcKadDd)5(Dc ), 

S(dn)  =  GAB8cGAB5(Dc)dn.  (56) 

Substituting  Eq.  51  and  Eq.  54- Eq.  56  into  the  left  side  of  Eq.  49  gives 

6  [  i,dn=  f  {PASA{dpa)  +  ZB[S(DA%B 
Jon  Jon 

+  [Qc  -  PAdB^dcN%  -  zB(gcJVA  -  5cKadDd)  (57) 
+  ( SAB  +  4>GAB)dcGAB]S(Dc)}dn. 

Two  applications  of  the  divergence  theorem  Eq.  19  and  repeated  integration  by  parts 
then  gives  the  following  equivalent  integral  form  of  Eq.  49: 

-  [  {[f„h  +  (PACCBC  +  SbPa)Dba]SV°  +  l Qc  -  Zg |B 

J'JJl 

-  ( ZaCEd  +  dBZA)DBA  -  ZB(dcNA  -  BcKadDd ) 

-  PA(dBdcVaDBA  +  dBpadcNB)  +  ( SAB  +  idGAB)dcGAB}5(Dc)}dn 

+  (f  [P*NA5va  +  ZBNB5(DA)}n=  (f  \pa5(fa  +  za5(Da)}Q. 

Jam  Jam 

(58) 

Assuming  this  global  equation  must  hold  for  admissible  variations  &c  and  5D,  local 
results  from  Eq.  58  are  the  Euler-Lagrange  equations  (i.e.,  force  balances)  and  the 
Neumann  boundary  conditions: 

dAPA  +  PbHab  -  PAHcbaFbA  +  PaNbCbc  +  (j PaCcbc  +  dBPA)dADB  =  0,  (59) 

d aZa  +  ZbHab  -  ZaHbc  +  8BZAdADB 

+  ZB(BCNA  -  BcKadDd  +  SaCbdDeb)  (60) 

+  PA(8B8c<paDBA  +  8Bipa8cNB)  -  ( SAB  +  i/jGab)8cGab  =  Qc; 

pa  =  PaNa,  za  =  ZbNb.  (61) 

Note  that  Eq.  47  has  been  used  in  Eq.  59  for  determining  the  (total)  horizontal 
covariant  derivative  of  the  2-point  tensor  P  =  PAdxa  <g)  t|j.  Equation  59  is  the 
(local)  balance  of  linear  momentum  for  quasi-statics.  Equation  60  will  be  referred 
to  as  the  (local)  balance  of  director  momentum  or  micromomentum. 

Balance  Eqs.  59  and  60  reduce  as  follows  when  configuration  spaces  are  (pseudo)- 
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Riemannian  ( G  independent  of  D): 


SaP‘  +  P^is  -  +  5BPA8ADB  =  0,  (62) 

8AZA  +  ZB~,AB  -  ZAjBc  +  8BZA8ADB  +  PB8Bdcv‘dADB  =  Qc ;  (63) 

(pseudo)-Minkowskian  (G  independent  of  A"  and  g  independent  of  x): 

aAP f  +  ( PbCcbc  +  8BPA)dADB  =  0,  (64) 

BaZb  +  8bZa8aDb  +  ZbCad8bDa  +  PA8B8c<p‘dADB 

(65) 

~(SAB  +  ijGAB)dcGAB  =  Qc; 

and  Cartesian  (global  metrics  Gab  =  At n  and  gab  =  8ab ): 

dAPA  +  8BPAdADB  =  0,  (66) 

dAZA  +  8BZAdADB  +  PA8B8c<padADB  =  Qc.  (67) 

The  model  framework  is  complete  upon  prescription  of  the  following  details.  For 
pseudo-Finsler  reference  space,  a  metric  tensor  G  is  introduced  over  the  domain  of 
interest  in  3,  from  which  all  connection  coefficients  are  derived  via  differentiation 
using  relations  listed  in  Section  2.1.  In  this  regard,  nonlinear  coefficients  can  be 
determined  from  the  spray  in  Eq.  9,  and  particular  forms  for  horizontal  and  vertical 
connection  coefficients  in  Eqs.  10  and  1 1  must  be  prescribed  (i.e.,  those  correspond¬ 
ing  the  Chern-Rund  connection  or  Cartan’s  connection).  The  main  requirement  for 
selection  of  horizontal  coefficients  is  that  (y/G)\A  must  vanish  for  the  form  of  di¬ 
vergence  theorem  in  Eq.  19  to  apply;  this  is  true  for  Chem-Rund  and  Cartan  con¬ 
nections.  Anaolgous  choices  must  be  prescribed  for  the  current  configuration  space, 
including  specification  of  metric  tensor  g ,  from  which  connection  coefficients  are 
derived  via  equations  in  Section  2.2.  Note  that  the  same  forms  of  metric  and  con¬ 
nections  need  not  be  prescribed  in  both  configurations;  that  is,  the  reference  space 
could  be  taken  as  Finslerian  and  the  current  configuration  space  Riemannian,  or 
vice-versa.  Regarding  deformation  kinematics  of  Section  2.3,  a  constitutive  equa¬ 
tion  may  be  added  for  specification  of  transformation  matrix  '0aR  in  Eq.  36,  but  this 
is  not  essential  for  the  solution  of  all  boundary  value  problems.  For  crystals,  a  con¬ 
venient  assumption  is  the  Cauchy-Born  rule,62-63  whereby  id\(X)  =  F'{  \X.  D(X)] 
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in  a  consistent  coordinate  basis.  A  particular  form  of  free  energy  function  tp  in  Eq. 
50  physically  appropriate  for  the  material  of  interest  must  also  be  invoked.  Then, 
given  prescribed  boundary  conditions  on  dDJl,  governing  Eqs.  59  and  60  represent, 
in  principle,  6  coupled  nonlinear  partial  differential  equations  for  6  unknown  fields 
ifa  \ X,  D(X)]  and  I)A{X).  If  Finsler  geometry  rather  than  pseudo-Finsler  geometry 
is  presumed,  then  a  fundamental  scalar  function  £  is  introduced  that  provides  G  by 
differentiation  via  Eq.  12  rather  than  direct  prescription,  and  similarly  for  1  and  g 
via  Eq.  31. 

A  few  major  differences  from  prior  literature  are  noted.  Specifically,  the  new  defi¬ 
nition  used  herein  for  the  deformation  gradient  (2-point)  tensor  in  Eq.  37 — a  delta- 
derivative — differs  from  that  in  reference  12 — a  partial  derivative — and  that  in  ref¬ 
erences  33  and  41 — a  covariant  derivative — and  can  be  interpreted  as  a  compromise 
between  the  other  latter;  this  compromise  further  enables  computation  of  Eq.  47. 
Certain  choices  or  options  for  metric  tensors  and  connection  coefficients  also  vary 
among  the  present  work  and  these  prior  works,  and  the  form  of  free  energy  function 
in  Eq.  50  differs  from  that  proposed  in  references  33  and  41. 

4.  Physics  of  Fracture 

The  first  of  3  boundary  value  problems  is  considered  in  the  present  section.  By 
construction,  these  problems  involve  fields  that  could  vary  in  only  1  or  2  rather 
than  all  3  spatial  directions.  Two  strategies  are  possible  for  formulating  the  pseudo- 
Finsler  theory  to  solve  such  problems:  either  define  a  general  3-D  theory  and  then 
reduce  the  equations  appropriately  for  the  problem  geometry,  or  construct  a  theory 
of  reduced  dimensionality  from  the  outset.  The  latter  approach  is  taken  in  Section 
4  to  obtain  physical  insight  while  maintaining  mathematical  simplicity. 

4.1  Problem  Geometry  and  Kinematics 

Considered  first  is  perhaps  the  simplest  physically  meaningful  application  of  the 
theory.  The  material  body  is  a  straight  1-D  bar  of  length  L0,  and  the  material  man¬ 
ifold  is  specified  as  (9DT  :  X  e  [0,  L0]}.  By  construction,  fields  vary  only  with 
X  =  X1  and  D  =  D1,  and  coordinates  A"2,  A"3, 1)1.  D3  are  superfluous.  A  Carte¬ 
sian  coordinate  system  suffices  for  {A"},  so  there  is  no  need  to  assign  a  metric 
tensor  with  dependence  on  X.  Consistent  with  these  protocols,  the  following  rela¬ 
tions,  which  are  reductions  of  more  general  definitions  and  identities  of  Section  2.1, 
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apply: 


X  =  X\  D  =  D1;  G  =  Gn(D),  Gn(D)  =  1  /G(D);  (68) 

7m  =  §0i  G  =  0,  7n  =  0;  G1  =  ^l^D1  =  0,  N ?  =  d.G1  =  0;  (69) 

Cm  =  \di  G  =  G'/2,  C\x  =  G'/{2G).  (70) 

The  reference  configuration  space  is  locally  Minkowskian.  Invoking  the  Chern- 
Rund  connection  with  vanishing  nonlinear  connection  coefficients  from  Eq.  69, 

Hi  =  Kl  =  rJi  =  ^qSiG  =  0;  Vl  =  Yl  =  0.  (71) 


For  the  current/deformed  configuration  of  the  bar,  with  deformed  material  manifold 
(m  :  x  G  |0.  L] }  ,  with  L  the  deformed  length  of  the  domain,  spatial  coordinates 
and  metric  components  are  of  the  assumed  Cartesian  form 

x  =  x1,  d  =  71;  g  =  g  n  =  1,  g11  =  l/g  =  1.  (72) 


All  spatial  connection  coefficients  of  Section  2.2 — linear  and  nonlinear — then  van¬ 
ish  identically. 

Motions,  deformations,  and  director  gradients  defined  in  Section  2.3  reduce  as  fol¬ 
lows  for  the  current  1-D  problem: 


x  =  <p(X,D),  d  =  7(X,  D);  D  =  D(X); 


(73) 


F(X,D)  =  F?(X,D)  =  5lX\X,D )  =  MgXD)  +  MqDD)  dDQX]  i  (74) 
J(X,D )  =  y/H, JG(D)F}(X,D )  =  [G(D)]~1/2[F(X,  D)];  (75) 

C(X,  D )  =  CU(X,  D )  =  F/(X,  D)guFl(X,  D )  =  [F(X,  F)]2;  (76) 

F|\  =  d\D  -  Nl  +  K]UD  =  dD/dX  =  D' .  (77) 


The  internal  state  variable  D  is  physically  identified  with  crack  opening  displace¬ 
ment.  Define  a  total  strain  measure  e  (which  includes  a  microdeformation  gradient 
contribution)  and  a  lattice  strain  measure  a  associated  with  stored  elastic  energy  as, 
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respectively, 


BD  B 

e  =  VC  -  1  =  F  -  1;  a  =  e-D'  =  e-  —  =  —(x  -  X  -  D ).  (78) 

Stretch  occurs  for  e  >  0,  contraction  for  e  <  0.  Also  introduce  a  constant  /  with 
dimensions  of  length  and  a  normalized  order  parameter  £: 

£  =  D/l,  £  =  D'/l;  a  =  e-  If.  (79) 

Constant  /  will  later  be  identified  as  the  value  of  crack  opening  displacement  at 
which  the  bar  supports  no  tensile  load.  Letting  k  denote  a  constant  depending  on 
the  material  of  interest,  a  more  specific  form  of  the  Minkowski  metric  in  Eq.  68  is 
introduced  as 

G(0  =  exp(2jfcf)  =>  G'/(2G)  =  Cln  =  k/l.  (80) 

The  length  of  a  referential  line  element  in  Eq.  14  and  the  corresponding  volume 
form  in  Eq.  15  become 

\dX\2  =  dX  ■  G  ■  dX  =  exp(2kD /l)dX  ■  dX,  dil  =  VGdX  =  exp(kD/l)dX, 

(81) 

such  that  expansion  occurs  when  k  >  0  and  contraction  when  k  <  0  iff;’  >  0.  The 
former  case  is  physically  representative  of  microscopic  dilatation  from  cracking  in 
crystalline  rocks,  for  example44,46  as  well  as  dilatation  from  core  fields  of  disloca¬ 
tions18,39,47  that  may  emerge  in  the  vicinity  of  crack  tips  in  more  ductile  crystals. 
Volume  changes  due  to  thermal  expansion  or  contraction  might  also  be  represented 
via  such  a  description.  Because  G  is  not  homogeneous  of  degree  zero  in  D,  the  ref¬ 
erence  configuration  space  is  not  strictly  of  Finsler  character,  but  rather  is  labeled  a 
pseudo-Finsler  space.16  Furthermore,  because  G  does  not  depend  explicitly  on  X, 
this  space  may  be  categorized  as  pseudo-Minkowskian. 

4.2  Energy,  Thermodynamic  Forces,  and  Balance  Laws 

Application  of  results  derived  in  Section  4.1  leads  to  the  following  forms  of  Eqs.  50 
and  52  for  the  present  1-D  example: 

^  =  i>(C,  D,  BD/BX,  G )  =  ip(F,  D,  BD/BX ,  G)  =  ^(e,  D,  BD/BX ,  G ).  (82) 
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Written  in  terms  of  lattice  strain  and  normalized  order  parameter,  this  becomes 

ip  =  iJ>[a(e,Z'),WMZ)].  (83) 

Motivated  again  by  phase  field  theory,43,60  the  following  sum  of  quadratic  forms  is 
invoked: 

V>  =  |A(l-£)V  +  T£2/(  +  T((£')2.  (84) 

Material  constants  are  defined  as  follows:  A  =  A  +  2/x  is  the  longitudinal  elastic 
stiffness,  T  is  the  crack  surface  energy  per  unit  reference  area,  and  regularization 
length  /  has  been  introduced  already  in  Eq.  79.  The  first  term  on  the  right  side  of  Eq. 
84  accounts  for  elastic  strain  energy  degraded  by  damage  associated  with  £  e  [0, 1], 
and  the  other  2  terms  combine  to  account  for  the  surface  energy  of  fracture.  This 
energy  function  contains  no  explicit  dependence  on  G. 


Nonzero  thermodynamic  forces  of  Section  3  are  then  obtained  by  direct  calculation 


Q  =  «‘  =  ^  =  7|  =  -f(1-«“2  +  2^  =  -^  +  2^  <86) 

(87) 

The  linear  momentum  balance  in  Eqs.  59  and  64  becomes 

9P(X,  D)  9P(X,  D)  3D  G’(D)  dD  AP  G>  AD 

gx  dD  ax  2G(D)  ax  iX  1G  AX  '  ’ 


The  micromomentum  balance  in  Eqs.  60  and  65  becomes 

dZ(X,D)  r  dZ(X,D)  G\D )  d2<p(X,D)  1  dD 

dX  l  dD  2 G{D)  dD 2  \  dX 


G\D) 

G(Dj 


1>(X,D)  =  Q{X,D). 


Substituting  from  Eqs.  80  and  85-87,  these  balance  laws  become 


^  =  -kP^- 
dX  dX' 


+  +2T/£"-2y£  =  2M^-T/(£')2]. 
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Relations  in  Eq.  90  constitute  a  pair  of  coupled  nonlinear  ordinary  differential  equa¬ 
tions  wherein  field  variables  P,  a,  £,  ip,  and  w  depend  ultimately  on  the  independent 
variable  X. 


4.3  Problem  Solutions:  Riemannian  Geometry 

Considered  first  is  the  simple  case  wherein  the  referential  configuration  space  is 
Riemannian  (in  fact,  a  Cartesian  structure),  with  k  —  0  — >  G  —  1  —  constant.  Even 
though  the  space  is  Cartesian,  the  kinematics  are  still  of  Finsler  character  since  p 
can  potentially  depend  on  both  D  and  X  rather  than  just  X  as  in  classical  continuum 
physics.  Balance  laws  in  Eq.  90  degenerate  to 


dP 

dX 


0; 


P 


a  _i_  d2yy 

'  r\  j-n  s 


1-e  ^ 


+  2T l£"  -  2j£  =  0. 


(91) 


The  first  of  Eq.  91  results  immediately  in  constant  stress  over  the  length  of  the  bar: 


P  =  P0  =  A(1  —  £)2a  =  constant. 


(92) 


The  solution  of  the  second  requires  further  specification  of  boundary  conditions. 
Two  particular  problems  corresponding  to  2  different  sets  of  boundary  conditions 
are  addressed  next:  homogeneous  damage  of  the  deformed  bar  over  [0,  L\  (i.e.,  mi¬ 
croscopic  fractures  evenly  distributed  along  the  length  of  the  bar)  and  localized 
damage  corresponding  to  a  globally  stress-free  deformed  state  (i.e.,  complete  ten¬ 
sile  fracture/rupture  of  the  bar). 

4.3.1  Homogeneous  Damage 

For  homogeneous  damage,  £'(A)  =  OVA"  e  [0,  L0\  =>-  £(0)  =  £(L0)  =  6 i- 
Boundary  conditions  on  displacement  are  prescribed  as  follows: 

<p{ 0,  D)  =  ip0  =  D  =  1£ht  ^(Lq,  D)  —  tpL  —  (1  +  a,H)L0  +  1£h-  (93) 

Here,  <pl  is  the  prescribed  coordinate  of  the  deformed  bar  at  x  =  L,  with  aH  and 
constants.  Equations  91  and  92  result  in 

P0  =  A(1  -  Zh  =  1/[1  +  2T/(Ma2H)\.  (94) 

Given  (pL ,  the  second  of  Eq.  93  and  Eq.  94  can  be  solved  simultaneously  for  the 
homogeneous  damage  field  £#,  stress  P0,  and  stretch  aH.  The  problem  kinematics 
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are  consistent  with  the  separable  decomposition  of  the  motion  (p  into 


<p[X,£(D)]  =  X(X)  +  l£(D),  F  =  dip/dX  =  y,  (95) 

where  for  homogeneous  damage  and  strain  fields, 

v[X,  £h]  =  (1  +  clh)X  +  l£H,  X  —  (1  +  uh)X,  F  —  1  +  aH.  (96) 
The  total  energy  of  the  degraded  elastic  bar  is  obtained  as 

'Hvl.(h)  =  >Ph  =  [  [A(1  -  ?„)24/2  +  reH/i]ix 

Jo  (97) 

=  [A(l  -  Z„)2a2H/2  +  T e„/l]L0. 

Remark:  Letting  tp(X,D)  — >  <p(X)  =  x(X)  and  modifying  Eq.  93  to  ip0  = 
0 ,(pL  =  (1  +  o,h)Lq  recovers  a  description  analogous  to  that  of  phase  field  the¬ 
ories  of  fracture  mechanics.43,64 

4.3.2  Stress-Free  State 

For  a  stress-free  state,  P  =  P0  =  OVA"  e  [0,  L0],  Boundary  conditions  on  the 
order  parameter  £(A")  or  internal  state  variable  representing  microdisplacement  are 
prescribed  as  follows: 


£(0)  =  D(Q)/l  =  1,  £(L0)  =  0.  (98) 

The  second  of  the  governing  equations  in  Eq.  91  becomes  the  homogeneous  linear 
second-order  ordinary  differential  equation  and  corresponding  exact  solution 

c  -  ai2  =  0  =s-  £(*)  =  eXp(y)/0  [exp(2.Y/0  -  exp(2io/0].  m 

1  —  exp(2L0/f) 

The  null  stress  condition  results  in  a(X)  =  0V£(A")  ^  1;  since  lattice  stretch  a( X) 
and  <p0  =  ip( 0,  D)  are  indeterminate  where  £  =  1,  the  total  displacement  at  the 
undamaged  end  of  the  bar,  (pi,  is  also  indeterminate.  Physically,  this  corresponds  to 
rigidly  displacing  the  bar,  once  disconnected  from  its  fully  fractured  site  at  X  =  0, 
without  altering  its  internal  energy.  Total  energy  is  the  integral 

*(0  =  =  r  m’)2l  +  £2/(]dX.  (100) 

7o 
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Remark:  Equation  99  contains  an  equation  and  associated  solution  identical  to  that 
attainable  through  phase  field  theories  of  fracture  mechanics.43-64 

4.4  Problem  Solutions:  Minkowskian  Geometry 

Considered  now  is  the  more  general  case  wherein  the  referential  configuration  space 
is  pseudo-Finslerian  (specifically,  a  locally  Minkowskian  structure)  with  k  ^  0  in 
Eq.  80.  Even  though  the  referential  metric  is  pseudo-Minkowskian,  deformation 
kinematics  are  again  of  Finsler  nature  since  x  =  tp(X,  D)  can  potentially  depend 
on  D  as  well  as  X .  Balance  laws  in  Eq.  90  apply;  the  first  results  in 

d P/P  =  -M£  =>P  =  P0  exp (— jfef),  (101) 

where  P0  is  the  constant  stress  corresponding  to  k  =  0  and/or  £  =  0.  The  solution  of 
the  second  balance  law  again  requires  further  specification  of  boundary  conditions. 
The  same  2  problems  considered  in  Sections  4.3.1  and  4.3.2  are  now  revisited  in  the 
context  of  the  pseudo-Minkowski  metric  in  Eq.  80,  recalling  that  k  >  0  accounts 
for  additional  microscopic  stretch  and  dilatation  in  the  damaged  zone  neglected  in 
the  Riemannian  metrical  representation  of  Section  4.3. 

4.4.1  Homogeneous  Damage 

For  homogeneous  damage,  £'(9f)  =  OVA"  e  [0,  L0\  =>-  £(0)  =  £(L0)  =  £#. 
Boundary  conditions  on  displacement  are  identical  to  Eq.  93: 

<p(0,  D)  =  (f0  =  D  =  /£# ,  <p{L0,  D)  =  ipL  =  (1  4-  o^)L0  +  1£h-  (102) 

Equations  90  and  101  result  in 

P  =  PH  =  Poexp(— k^H)  =  A(1  —  £fl-)2Off  =  constant, 

Aa^(l  -  £#)  -  2(Y /1)£h  =  Skip  (an,  £h)- 

Given  tpL,  Eqs.  102  and  103  can  be  solved  simultaneously  for  the  homogeneous 
damage  field  £H,  stress  P0,  and  stretch  aH.  The  problem  kinematics  are  consistent 
with  the  separable  decomposition  of  Eqs.  95  and  96.  The  total  energy  of  the  elastic 
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bar  is  obtained  as 


n<PL,tH)  =  ^h=  A  A(1  -  iufa\l2  +  T&/J]dX 

Jo  (104) 

=  [A(l-^)2a^/2  +  T^//]L0, 

identical  in  form  to  Eq.  97  but  with  potentially  different  values  of  £#  and  a  rr  for  a 
prescribed  ip l  when  k  is  nonzero. 

Shown  in  Fig.  1  are  £  =  £#,  P  =  PH,  and  T  =  versus  applied  tensile  displace¬ 
ment,  computed  via  Eqs.  103  and  104.  Representative  material  parameters  are  taken 
as  A  —  n  —  A/3  =  10nN/m2,  T  =  lN/m,  and  l  =  lCT9m.43-60  The  domain  size  is 
L0  =  103/,  and  the  Weyl  scaling  parameter  k  is  varied  from  0  to  In  2,  with  the  latter 
(maximum  considered)  value  corresponding  to  a  maximum  volume  form  scaling  of 
y/G  =  2  at  £//  =  1.  Since  resulting  volume  changes  are  considerably  large,  the 
Weyl  scaling  could  be  interpreted  as  giving  rise  to  a  fictitious  damaged  configura¬ 
tion,  similar  to  that  envisioned  in  nonlinear  continuum  damage  mechanics.65  For 
any  fixed  value  of  k,  the  order  parameter  £  increases  monotonically  with  increasing 
tensile  displacement  (Fig.  la),  tensile  stress  P  increases  to  a  maximum  and  then  de¬ 
creases  (Fig.  lb),  and  energy  T  increases  monotonically  (Fig.  lc).  As  k  increases, 
the  value  of  £  tends  to  decrease  for  0  <  (<Pl/Lq  —  1)  <  0.5,  while  P  and  T  tend  to 
increase  over  the  same  range  of  displacement.  Notably,  as  evident  from  Fig.  lb,  the 
peak  stress  and  the  applied  displacement  at  which  the  peak  stress  is  reached  both 
increase  significantly  with  increasing  k,  implying  an  increase  in  tensile  strength  and 
stability  of  the  material  commensurate  with  microscopic  dilatation  represented  by 
positive  values  of  k.  The  increases  in  stress  and  energy  correlate  with  a  decrease 
in  order  parameter  since  both  P  and  strain  energy  density  contain  a  multiplication 
factor  of  (1  —  £)2.  The  contribution  of  the  /£#  term  in  Eq.  96  is  negligible  since 
I/Lq  <C  1;  therefore,  differences  between  the  present  Finsler  solution  with  k  —  0 
and  phase  field  theory43-64  are  also  negligible. 
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(a)  (b) 


(c) 

Fig.  1  Tensile  deformation,  homogeneous-state  solutions,  l/L0  =  10-3:  (a)  order  parameter 
£  =  D /l,  (b)  normalized  tensile  stress,  and  (c)  normalized  total  energy 


4.4.2  Stress-Free  State 

As  in  Section  4.3.2,  for  a  stress-free  state  P  =  OVA"  e  [0,  L0\.  Boundary  conditions 
on  the  order  parameter  £( X)  are  prescribed  as  in  Eq.  98: 

f  (0)  =  D(0)/l  =  1,  f  (L0)  =  0.  (105) 

The  second  of  governing  equations  in  Eq.  90,  with  Eq.  84,  becomes  the  nonlinear 
second-order  ordinary  differential  equation 

f  -  ?/i2  +  fc[(€')2  -  !/>/( TO]  =  o  f"  =  (?/i2)(l  +  kO-  (106) 
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Defining  £  =  £'  such  that  £"  =  £  •  d£/d£,  this  can  be  transformed  into  the  nonho- 
mogeneous  first-order  differential  equation  and  corresponding  general  solution 

CdC  =  (£//2)(l  +  fcO d£  =>  c  =  ±(e//)Vl  +  2A:e2/3  +  c1/e.  (107) 

In  this  case,  integration  constant  c\  vanishes  and  the  negative  root  applies.  The 
second  of  Eq.  107  can  be  integrated  to  give  the  following  implicit  solution  for  £(X), 
which  is  then  evaluated  via  numerical  quadrature: 

d£  =  -K/i)^l  +  2*^/3 d.Y  =>  X(()  =  f  ~id“  (108) 

J i  «yl  +  2/ccr/3 

The  null  stress  condition  again  results  in  a(X)  =  0V£(X)  ^  1  so  that  <p0  and 
remain  indeterminate.  Total  energy  is  the  integral 

*(o  =  =  lU)  rn'Yi + e/i]  dx.  (io9) 

Jo 

The  total  energy  per  unit  cross-sectional  area  of  the  bar  computed  via  Eq.  109  is 
shown  in  column  2  of  Table  1,  where  L0  =  1  for  normalization.  This  energy  T/r 
increases  slightly  with  increasing  k,  with  the  k  —  0  solution  identical  to  the  stress- 
free  solution  from  phase  field  theory.43,64  A  value  of  =  T  corresponds  to 

Griffith’s  theory  of  mode  I  brittle  fracture,  recalling  that  material  property  T  is 
surface  energy.  Shown  in  Fig.  2  are  profiles  of  £  computed  via  Eq.  108,  with  a 
domain  size  of  L0  =  10 1  and  the  same  range  of  Weyl  scaling  factor  k  considered  in 
Section  4.4.1.  Regardless  of  k,  the  value  of  £  drops  off  rapidly  from  its  maximum  at 
X  =  0  with  increasing  X.  Increasing  k  results  in  a  small  decrease  in  £  for  X  <  L0, 
and  an  increase  in  |£'|  (i.e.,  a  sharper  fracture  profile). 

Table  1  Stress-free  1-D  solutions  for  l/L0  =  0.1:  total  energy 


Weyl  Scaling  Factor 

Tension/Shear:  ^!f/(TLq) 

k  =  0 

1.0091 

k  =  In  | 

1.0100 

k  =  In  | 

1.0117 

k  =  In  | 

1.0138 

k  =  In  2 

1.0159 
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Fig.  2  Axial  or  shear  stress-free  solutions,  l/L0  =  0.1:  (a)  0  <  X  <  L0  and  (b)  £:  0  <  X  < 

O.ILq 


5.  Physics  of  Slip 

The  second  of  3  boundary  value  problems  is  considered  in  this  section,  involving 
simple  shearing  of  a  2-D  nonlinear  elastic  slab.  Herein,  a  general  free  energy  func¬ 
tion  is  postulated,  applicable  for  any  3-D  deformation  modes,  and  then  specified  to 
the  present  geometry  and  kinematics. 

5.1  Problem  Geometry  and  Kinematics 

The  material  body  is  an  elastic  slab  of  length  L0  and  infinite  width  and  thickness.  In 

2  dimensions,  the  material  manifold  is  specified  as  {921 :  X 1  e  [0,  L0\,  |X2|  e  oo}. 
Regarding  the  third  (out-of-plane)  direction,  plane  strain  conditions  are  imposed. 
The  internal  state  vector  is  restricted  as  {DA}  — >  {0,  D2 ,  0}.  By  construction,  fields 
vary  only  with  A"  =  X1  and  D  =  D2,  and  coordinates  A"3,/}1,  D3  are  superflu¬ 
ous.  A  Cartesian  coordinate  system  suffices  for  {AT}  so  metric  tensor  G  contains 
no  dependence  on  A".  Consistent  with  these  protocols,  the  following  reductions  of 
definitions  and  identities  of  Section  2.1  hold: 


{X,F}  =  {X\X2},  D  =  D2- 


(HO) 


G(D) 


Gu(D)  G  12(D) 
G12(D)  G22(D) 


G 


G11G22 


G 


2 

12) 


(HD 
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7 abc  —  ^{BaGbc  +  BbCac  —  9cGab)  —  0,  GA  —  \^CDBDC  —  0;  (112) 
NA  =  dBCrA  =  0  =>  SA(-)  =  dA(-y,  (113) 

<3*111  =  <3*122  =  <3*212  =  0,  C'222  =  G'2-2/2,  <3*112  =  ~G'u/2, 

(114) 

<3*121  =  <3*211  =  <3*ii/2,  <3*221  =  <3*12- 

The  reference  configuration  space  is  locally  Minkowskian.  Invoking  the  Chern- 
Rund  connection  with  vanishing  nonlinear  connection  coefficients  from  Eq.  112, 


HAC  =  Kin  =  r 


BC 


BC 


=  \G 


AD 


( ScGbd  +  5BGcd  —  $dGbc)  —  0; 


VA  —  vA  —  n 

V  BC  —  1  un  —  u 


(115) 


BC 


For  the  current/deformed  configuration  of  the  slab,  with  deformed  material  mani¬ 
fold  (m  :  x 1  G  [0,  L\,  \x2\  G  oo},  where  L  is  the  deformed  length  of  the  domain, 
spatial  coordinates  and  metric  components  are  of  the  Cartesian  form 

{x,y}  =  {xl,x2},  d  =  d2]  gab  =  Sab,  9  =  1-  (116) 


All  spatial  connection  coefficients  of  Section  2.2 — linear  and  nonlinear — then  van¬ 
ish  identically. 

Motions,  deformations,  and  director  gradients  defined  in  Section  2.3  reduce  as  fol¬ 
lows  under  simple  shear,  with  (p  =  v  +  Y  and  e  denoting  deformation  and  strain  in 
the  shearing  (Y)  direction: 


x  =  X,  y  =  p(X,  Y,  D )  =  Y+v(X,  D);  d  =  0(X,  D);  D  =  D{X);  (117) 

dx(X)  dx(X) 


F(X,  D)  = 


dX 


dY 


dip(X,Y,D)  .  dip(X,Y,D)  dD(X)  d^p{X,Y,D) 

dX  '  dX  dX  dY 


d  v[X,D(X)\ 
dX 


ox 

OX 

i  o' 

e(X,D)  1 

J(X,  D)  =  ys/G(D)Fi(X .  D)FpX,  D)  =  \G(D)]-'!2  =  J(D); 
<3*ii  =  1  +  e2,  <3*12  =  <3*2i  =  <3*22  =  1;  det  C  —  1; 

Df,  =  f)J)  -  N2  +  K22D  =  8D/8X  =  I)’. 


(US) 

(119) 

(120) 
(121) 
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The  internal  state  variable  D  is  physically  identified  with  a  slip  discontinuity  or 
slipped  displacement.  The  defect  associated  with  D  could  be  a  shear  band  in  a 
(poly)crystalline  metal,66  a  stacking  fault  in  a  crystal  lattice,19  or  a  mode  II  crack  in 
a  brittle  solid.43  It  becomes  useful  to  define  a  lattice  shear  strain  measure  a ,  defor¬ 
mation  gradient  F,  deformation  tensor  C,  and  Jacobian  determinant  J  as  follows: 


a  =  e  ±  D' 


d 

dX 


(: y-Y±D ), 


F(X,D,D') 


1  0 
a(X,  £>,£>')  1 


(122) 


C\  i  —  1  +  a2,  C  12  —  C 2  i  —  a,  C22  —  1;  J  —  \J  det(C,As)-  (123) 

As  in  Section  4,  introduced  are  a  constant  l  with  dimensions  of  length  and  a  nor¬ 
malized  order  parameter  £  e  [0, 1]: 


£  =  ZJ//,  £'  =  £>'//;  a  =  e±Z£'.  (124) 


Constant  l  will  be  identified  as  the  value  of  shear  slip-displacement  at  which  the  slab 
supports  no  shear  stress.  Letting  k  denote  a  constant  depending  on  the  material,  a 
more  specific  form  of  the  Minkowski  metric  in  Eq.  1 10  is  invoked: 


G{D) 


VG(D)  0 
0  VG(D)  ’ 


G(0 


exp  (2  fc£)  => 


G 

2G 


G\ 


ii 


G 


ii 


G22 


_  k 
~  ~l 
(125) 


This  can  be  viewed  as  a  Weyl  transformation  or  Weyl  rescaling48  of  the  Cartesian 
metric  5ab-  Also  used  later  is  the  second  component  of  the  trace  of  Cartan’s  tensor 
ofEq.  114: 


CfA  =  GabC-2ab  =  l(GuG'u  +  G22G’22  +  4  G12G'12)  =  k/l.  (126) 


The  length  of  a  referential  line  element  in  Eq.  14  and  the  corresponding  volume 
form  in  Eq.  15  become 


IdXl 


(dX  •  dX  +  d Y  ■  dF), 


dVt  =  exp 


dXAdX.  (127) 


For  £  >  0,  expansion  occurs  when  k  >  0  and  contraction  when  k  <  0.  Physical 
justification  follows  similar  arguments  as  given  in  Section  4.1:  shear  fractures  may 
result  in  dilatation  as  rough  crack  faces  slide  over  one  another,67  while  (full)  dis- 
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locations  or  partial  dislocations  associated  with  shear  bands  or  stacking  faults  may 
result  in  local  dilatation  due  to  nonlinear  elastic  and  core  effects. 68-70  Because  G 
is  not  homogeneous  of  degree  zero  in  D,  the  reference  configuration  space  is  not 
strictly  of  Finsler  character,  but  is  a  pseudo-Finsler  space,16  and  because  G  does  not 
depend  on  X,  this  space  may  be  further  categorized  as  pseudo-Minkowskian. 

5.2  Energy,  Thermodynamic  Forces,  and  Balance  Laws 

For  a  compressible  neo-Hookean  elastic  solid  with  strain  energy  function  W  de¬ 
pending  on  lattice  deformation  tensor  C,  the  following  general  form  of  total  free 
energy  density  in  Eq.  52  is  assumed: 

V>( CAB ,  DA,  DA,  GAB )  =  W[CAB(CAB ,  DA),  DA]  +  f[DA,  DA,  Gab],  (128) 

where  function  /  accounts  for  surface  energy.  Specifically,  extending  prior  phase 
field  theory43  to  Finsler-geometric  continuum  mechanics,  let  /i  and  A  denote  the 
usual  isotropic  elastic  constants  and 

W  =  [\ii{Cab5ab  ^\nJ +\\{\nJ)2]{l-^)\  (129) 

/  =  (T/l)(DADB5AB/l2  +  DAB8AC5BDDfD).  (130) 

A  stress  tensor  P  associated  with  W  is 

*?~Wx~Wx~  W**6*6**  +  (AlnJ  -  "FAKi  -  O’-  (131) 

For  the  present  case  of  simple  shear  with  Eqs.  120  and  123  now  applied, 

D')}}/dFbB  =  (da/at)%6%  =  6"b6%  ^  P, ?  =  P ?.  (132) 

This  identity,  with  Ff  =  —  (  F- 1  =  a  and  J  —  1  in  Eq.  131,  results  in  the  only 

nonzero  stress  components 

P  =  P21  =  P12  =  M  l-02a.  (133) 

The  balance  of  angular  momentum  in  Eq.  53  is  verified  as  satisfied  since,  for  the 
present  problem, 

aob  =  crba  ^  Pi  =  P\  -  eP\  =  P\.  (134) 
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Derivations  from  Section  5.1  result  in  the  following  form  of  Eq.  128  that  applies  for 
the  present  example: 


ip  =  ip(e,  D,  dD/dX)  =  iP( dv/dX ,  D,  dD/dX) 

=  W(e,  D,  dD/dX)  +  f(D,  dD/dX). 

Written  in  terms  of  lattice  shear  strain  and  normalized  order  parameter,  this  be¬ 
comes 

^  =  iP[a(e,0,p,p'}  =  W[a(e,0^}  +  /(£,£')■  (136) 

Adding  /  to  the  elastic  strain  energy  W  reduced  under  simple  shear  leads  to 

^  1  -  0  v  +  t  e/1  +  ri(02-  (137) 


The  first  term  on  the  right  side  of  Eq.  137  accounts  for  elastic  strain  energy  degraded 
by  damage  associated  with  £  e  [0,1],  and  the  other  2  terms  combine  to  account 
for  surface  energy  associated  with  the  particular  class  of  shear  defect  (shear  band, 
stacking  fault,  mode  II  crack,  etc.)  under  consideration.  This  energy  function  con¬ 
tains  no  explicit  dependence  on  G  and  is  nearly  identical  to  that  of  Eq.  84,  differing 
only  in  the  elastic  constant  (//  versus  A  =  A  +  2/i)  and  physical  meanings  of  the 
lattice  strain  variable  a,  the  order  parameter  £,  and  possible  values  of  T  and  /. 


Relevant  (i.e.,  possibly  nonzero)  thermodynamic  forces  of  Section  3  are  then  ob¬ 
tained  by  direct  calculation  as 


P  =  Pl 


dip 

Ik 


dW  da(e,  £') 
da  de 


M1  -02°; 


(138) 


n-n  _  dip  _ldip  _  ^  2  ,  J 

Q  dD  i  dp  i^1  +2i2^ 


■W=3+2h  (139) 


(140) 


Notice  that  Eq.  138  is  consistent  with  Eq.  133.  The  linear  momentum  balance  in 
Eqs.  59  and  64  becomes 


dP(X,  D )  dP(X,  D)  dD  „  G'(D)  dD 
dX  +  dD  dX  +P2G(D)dX 


d P  G'  dD 
dX+P2GdX 


0.  (141) 
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Micromomentum  balances  in  Eqs.  60  and  65  become 


d Z(X,  D) 
dX 


dZ(X,D)  +zG'(D) 


dD 


2  G(D) 


+  P 


d2v(X,  D ) 


dD 2 
G\D) 


dD 

dX 


G(D) 


i/j(X,D)  =  Q(X,D). 


(142) 


Substituting  from  Eqs.  125,  126,  and  138-140,  these  balance  laws  become,  respec¬ 
tively, 


^  =  -kP^- 
dX  dX' 


P 


a  d2v 
H - £ 

1  o  >-o  S 


i-e  de 


+  2Yl£"-2j£  =  2k[^-Yl(£f)2}. 


(143) 

The  relations  in  Eq.  143  are  2  coupled  ordinary  nonlinear  differential  equations 
wherein  field  variables  P,  a,  £,  v  =  p  —  Y,  and  ru>  depend  ultimately  on  independent 
variable  X. 


Remark:  Balance  laws  in  Eqs.  90  and  143  are  mathematically  identical,  but  phys¬ 
ically  represent  2  different  problems  (i.e.,  tensile  deformation  and  simple  shear  de¬ 
formation),  with  P  representing  axial  stress  in  the  former  and  shear  stress  in  the 
latter. 


5.3  Problem  Solutions:  Riemannian  Geometry 

Considered  first  is  a  Riemannian  referential  configuration  space  (in  fact,  a  Cartesian 
structure),  with  k  —  0  — >  G  —  1  —  constant.  As  in  Section  4.3,  kinematics  are  still 
Finslerlian  since  99  can  potentially  depend  explicitly  on  both  D  and  X  rather  than 
just  A".  Balance  laws  in  Eq.  143  reduce  to 


dP 

dX 


0; 


P 


a 


+ 


dh^; 

de\ 


+  2YlC-2j£ 


0. 


(144) 


The  first  of  Eq.  144  results  immediately  in  spatially  constant  shear  stress  compo¬ 
nents: 

P  =  P0  =  /x(l  —  £)2a  =  constant.  (145) 

The  solution  of  the  second  balance  equation  requires  more  precise  boundary  con¬ 
ditions.  Two  problems  corresponding  to  2  different  sets  of  boundary  conditions  are 
addressed:  homogeneous  damage  of  the  deformed  slab  over  A"  e  [0,  L\  (i.e.,  mi¬ 
croscopic  shear  fractures  or  slip  bands  D  =  D2  evenly  distributed  along  the  finite 
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length  of  the  slab)  and  localized  damage  corresponding  to  a  globally  stress-free 
deformed  state  (i.e.,  complete  shear  failure  of  the  slab  along  a  plane  X  =  constant). 

5.3.1  Homogeneous  Damage 

For  homogeneous  damage,  £'(9f)  =  OVA"  e  [0,  L0\  =>  £(0)  =  £(L0)  = 
Boundary  conditions  on  shear  displacement  (or  on  p  along  Y  =  0)  are  prescribed 
as  follows: 


u(0,  D)  =  Vo  =  ±D  =  ±l£H,  v(L0,  D )  =  vL  =  aHLo  ±  1£h-  (146) 

Here,  v jJ  is  the  prescribed  displacement  of  the  deformed  slab  at  x  =  L,  with  an  and 
£#  constants.  Evolving  rigid  body  displacement  (i.e.,  shear  slip)  at  the  end  of  the 
slab  is  quantified  by  vq.  Equations  144  and  145  result  in 

P0  =  //(l  -  6r)  V,  Zh  =  1/[1  +  2T /(iila2H)\.  (147) 

Given  vl,  Eqs.  146  and  147  can  be  solved  simultaneously  for  the  homogeneous 
damage  field  £H,  stress  P0,  and  lattice  shear  strain  an-  The  problem  kinematics  are 
consistent  with  the  separable  decomposition  of  the  motion  p  into 

<p[X,  £(D),  Y]  =  v(X,  D )  +  Y  =  X(X)  +  Y  ±  ^(D),  F *  =  dp/dX  =  X' , 

(148) 

where  for  homogeneous  damage  and  strain  fields, 

<p[X,ZH,Y]  =  aHX  +  Y±ltH,  X  =  aH  X,  F?  =  aH.  (149) 

The  total  energy  of  the  slab  per  unit  width  in  the  V'-dircction  is  obtained  as 

*(vl,Zh)  =  *h  =  /  V(  1  -  infap 2  +  T^/J]dX 

Jo  (150) 

=  Ml  -  (h)Vh/ 2  +  T(l ,/l]La. 

Remark:  Letting  p(X,  Y,  D)  — >  p(X,  Y)  =  x{X)  +  Y  and  modifying  Eq.  146 
to  =  0,  vl  =  a  n  Lo  recovers  a  description  identical  to  a  prior  geometrically 
nonlinear  phase  field  study.43 
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5.3.2  Stress-Free  State 


For  a  stress-free  state,  P  =  P0  =  OVA"  e  [0,  L0].  Boundary  conditions  on  the  order 
parameter  £(A")  or  internal  state  variable  representing  microshearing  are  prescribed 
as  follows: 

m  =  D{0)/l  =  ll  £{L0)  =  0.  (151) 

The  second  of  governing  equations  in  Eq.  144  becomes 

£"  -  £//2  =  0  =>  £(AT)  =  ci  exp (X/l)  +  c2  exp (-X/1).  (152) 

Constants  c\ ,  c2  are  determined  by  boundary  conditions  in  Eq.  151,  leading  to  a 
particular  analytical  solution  identical  to  that  in  the  second  of  Eq.  99.  The  null 
stress  condition  results  in  a( X)  =  0V£(A")  ^  1;  since  lattice  shear  strain  a( X) 
and  shear  displacement  v0  =  v(0,  D)  are  indeterminate  where  £  =  1,  the  total 
shear  displacement  at  the  undamaged  end  of  the  slab,  vl,  is  also  indeterminate. 
Physically,  this  corresponds  to  rigidly  displacing  the  bar,  once  disconnected  from 
its  fully  localized  and  degraded  site  at  X  =  0,  without  altering  its  internal  energy. 
Total  energy  per  unit  width  is  the  integral 

* (o  =  *F=  [ 0  m')2i + e/i]dx.  (153) 

Jo 

Remark:  Eq.  152  is  a  result  identical  to  that  of  a  phase  field  theory  of  shear  failure43 


5.4  Problem  Solutions:  Minkowskian  Geometry 

Considered  now  is  a  referential  configuration  space  that  is  pseudo-Finslerian  with 
k  ^  0  in  Eq.  125.  Even  though  the  referential  metric  is  pseudo-Minkowskian,  de¬ 
formation  kinematics  are  again  of  Finsler  nature.  Balance  laws  in  Eq.  143  apply; 
the  first  results  in 

d P/P  =  -kd£  =>P  =  P0  exp (— fcf),  (154) 

where  P0  is  the  constant  shear  stress  corresponding  to  k  —  0  and/or  £  =  0.  The 
2  problems  considered  in  Sections  5.3.1  and  5.3.2  are  revisited  in  the  context  of 
the  pseudo-Minkowski  metric  in  Eq.  125,  recalling  that  k  >  0  accounts  for  dilata¬ 
tion  in  the  damaged  or  intensely  sheared  zone  omitted  in  the  Riemannian  metrical 
representation  of  Section  5.3. 
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5.4.1  Homogeneous  Damage 

For  homogeneous  damage,  £'( X)  =  OVX  e  [0,  L0\  =>•  £(0)  =  £(L0)  =  £/?. 
Boundary  conditions  on  displacement  v  =  if  —  Y  are  identical  to  Eq.  146: 


f(0 ,D)  =v0  =  ±D  =  ±l£H,  v(L0,  D)  =  vL  =  aHL0  ±  l£H.  (155) 


Equations  143  and  154  give 

P  =  P0exp(— k£H)  =  A*(l  —  Gr)2°ir  =  constant, 
/t^(  1  —  6r)  —  2(T /1)£h  =  2ki/j(aH,  £#). 


(156) 


Given  vL,  Eqs.  155  and  156  can  be  solved  simultaneously  for  £H,  P0,  and  aH- 
Kinematics  are  consistent  with  separable  decompositions  in  Eqs.  148  and  149.  The 
total  energy  per  unit  width  of  the  elastic  slab  is 


* K6r)  =  =  [L°H  1  -  t,H?a2H/ 2  +  TZ2H/l]dX 

Jo 

=  [Ml-^)2aV2  +  T^//]L0. 


(157) 


Shown  in  Fig.  3  are  ^  =  £h,  P  =  Ph,  and  T  =  T h  versus  applied  shear  displace¬ 
ment,  computed  via  Eqs.  156  and  157,  where  the  positive  choice  ip(X,Y,D )  = 
x(X)  +Y  +  D  is  applied  in  Eqs.  148  and  155.  Representative  material  parame¬ 
ters  are  identical  to  those  invoked  in  Section  4.4.1:  /i  =  109N/m2,  T  =  lN/m, 
l  =  10_9m,  L0  =  103/,  and  0  <  k  <  In  2.  Trends  are  similar  to  those  depicted  in 
Fig.  1.  For  fixed  k,  £  increases  monotonically  with  increasing  shear  displacement 
(Fig.  3a),  shear  stress  P  increases  to  a  maximum  and  then  decreases  (Fig.  3b),  and 
energy  T'  increases  monotonically  (Fig.  3c).  As  k  increases,  £  tends  to  decrease  for 
0  <  vl/Lq  <  0.5,  while  P  and  T  tend  to  increase.  Effects  of  k  are  more  pronounced 
in  Fig.  3  than  in  Fig.  1,  with  the  difference  due  to  different  elastic  constants  for  shear 
and  uniaxial  tension  (i.e.,  //  and  A  =  3/r).  Peak  shear  stress  and  applied  displace¬ 
ment  at  which  peak  stress  is  attained  both  increase  significantly  with  increasing  k, 
implying  an  increase  in  shear  strength  and  stability  of  the  material  commensurate 
with  microscopic  dilatation  represented  by  k  >  0.  Again,  increases  in  stress  and 
energy  correlate  with  a  decrease  in  order  parameter  since  both  P  and  strain  energy 
density  W  are  affected  by  a  multiplication  factor  of  (1  —  £)2:  see  Eqs.  129  and  133. 
The  contribution  of  the  Finsler  kinematic  /£#  term  in  Eq.  149 — in  other  words,  rigid 
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slip — is  negligible  for  l/L0  -C  1,  so  the  present  Finsler  solution  with  k  —  0  and 
phase  field  theory43  provide  nearly  identical  results  when  microscopic  dilatation  is 
omitted. 


(c) 

Fig.  3  Shear  deformation,  homogeneous-state  solutions,  l/L0  =  10-3:  (a)  order  parameter 
£  =  D /l,  (b)  normalized  shear  stress,  and  (c)  normalized  total  energy 


5.4.2  Stress-Free  State 

Just  as  imposed  in  Section  5.3.2,  let  P  =  OVA"  e  [0,  L0\.  Boundary  conditions  on 
£(A)  are  prescribed  as  in  Eq.  151: 


£(0)  =  D(0)/l  =  1,  £(£o)=  0.  (158) 
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The  second  of  governing  equations  in  Eq.  143,  with  the  reduced  form  of  energy 
density  in  Eq.  137,  becomes  the  nonlinear  second-order  ordinary  differential  equa¬ 
tion 

e  -  ai2  +  m')2  -  V’/to]  =  o=>t"  =  (^//2)(  i + k$.  (159) 

Defining  (  =  £'  such  that  —  £  ■  d£/d£,  this  can  be  transformed  into  the  nonho- 
mogeneous  first-order  differential  equation  and  corresponding  general  solution 

cdc  =  (£//2)(i  +  ho#  =>  c  =  ±(z/iwi  +  2ke/3  +  Cl/e.  ci60) 

As  in  Section  4.4.2,  the  latter  can  be  rewritten  and  then  integrated  to  give 

d?  =  -mV TTWTM  -  a-K)  =  /  (161) 

The  null  stress  condition  again  results  in  a(X)  =  0V£(X)  ^  1  with  v0  and  vL 
indeterminate  due  to  admissible  rigid  body  motion.  Total  energy  per  unit  width  is 
the  line  integral 

* (o  =  *f=  [ 0  mn + e/i] dx.  (i62) 

7o 

Remark:  Developments  in  Section  5.4.2  are  mathematically  identical  to  those  of 
Section  4.4.2,  but  with  different  physical  implications  (i.e.,  shear  failure  versus  ten¬ 
sile  failure).  Solutions  depicted  in  Fig.  2  and  Table  1  apply  here,  with  identical 
mathematical  trends  to  those  discussed  in  Section  4.4.2,  including  close  agreement 
of  shear  failure  or  mode  II  fracture  energies  among  the  present  theory,  phase  field 
theory,  and  Griffith’s  theory.  Physically,  the  pseudo-Finsler  theory  predicts  that  di¬ 
latation  is  associated  with  an  increase  in  slip  strength  or  crack  sliding  resistance 
as  well  as  failure  energy.  Such  predictions  agree  with  physical  observations.67  In 
an  intensely  sheared  zone,  effects  of  microscopic  friction  or  locking  of  asperities 
increase  when  the  material  dilates  in  such  a  zone.  Recall  from  Section  5.1  that 
the  dilatation  here  may  be  the  result  of  heterogeneous  (i.e.,  imperfect)  microscopic 
fractures  (as  opposed  to  perfect  cleavage),  expansion  of  the  material  due  to  nonlin¬ 
ear  elastic  effects  associated  with  defect  cores  (e.g.,  dislocation  cores,  point  defects, 
and  stacking  faults68),  and/or  thermal  expansion  due  to  temperature  rise  in  adiabatic 
shear.66  It  should  be  noted,  however,  that  the  present  variational  model  does  not  ex¬ 
plicitly  monitor  time-dependent  dissipated  energy  associated  with  sliding  friction 
or  dislocation  glide. 
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6.  Physics  of  Cavitation 


The  final  of  3  boundary  value  problems  involves  radial  expansion  of  a  spherical 
nonlinear  elastic  body.  Cavitation  (i.e.,  void  or  vacancy  formation  and  growth)  may 
occur  uniformly  throughout  the  domain  and/or  localized  at  or  very  near  to  the  cen¬ 
ter  of  the  sphere,  depending  on  particular  boundary  conditions  imposed  in  what 
follows.  A  general  free  energy  function  is  postulated  as  in  Section  5  and  then  spe¬ 
cialized  to  the  current  geometry,  but  here  curvilinear  (specifically,  spherical)  coordi¬ 
nates  are  needed.  This  complicates  the  present  analysis  relative  to  those  considered 
in  Sections  4  and  5  that  were  tractable  via  Cartesian  frames.  The  metric  tensor  nec¬ 
essarily  depends  on  position  and  may  also  depend  on  the  internal  state  vector;  when 
the  latter  applies,  a  completely  pseudo-Finslerian  referential  configuration  space 
results. 

6.1  Problem  Geometry  and  Kinematics 

The  material  body  is  an  elastic  sphere  of  radius  R0.  The  referential  material  mani¬ 
fold  is  specified  as  (SOT  :  R  =  X1  G  [0,  R.q]  ,  0  =  X2  E  [0, 7r],  $  =  X3  E  (— n,  7r]}. 
By  construction,  spherical  symmetry  conditions  are  imposed  so  that  solution  field 
variables  do  not  depend  on  angular  coordinates  0,  <f>.  The  internal  state  vector  is 
restricted  to  have  a  radial  component  only:  {DA}  — >■  { D 1,  0,  0}.  Metric  tensor  G 
necessarily  depends  on  X  and  possibly  depends  on  D.  Consistent  with  these  proto¬ 
cols,  definitions  and  identities  of  Section  2.1  result  in 

{R,  0,  =  {X1,  X2,  X3},  {D\  D2,  D3}  =  {D,  0,  0}.  (163) 

Denoting  by  G  =  G(X)  the  usual  metric  tensor  of  Euclidean  space  in  spherical 
coordinates25-71  and  B  =  B ( I) )  a  differentiable  scalar  function  of  the  internal  state, 
the  following  separable  form  of  the  metric  tensor  G(X,  D )  is  assumed: 

G{X,D)  =  G{X)B(D ) 


'Gu 

0 

0 

'l 

0 

0 

0 

G22(X) 

0 

B{D)  = 

0 

R2 

0 

_  0 

0 

G33(X)_ 

0 

0 

R2  sin2  0_ 

(164) 

G  =  B3G  =  B3R 4  sin2  0.  (165) 
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Let  denote  Christoffel  symbols  of  the  second  kind  derived  from  G,  and  let  7  BC 
denote  those  derived  from  G.  From  Eq.  6,  these  are  equivalent: 

7  bc  =  Gad^bcd  =  (B  1Gad)(B^bcd)  =  7  bc-  (166) 

Nonzero  Christoffel  symbols  resulting  from  Eq.  164  are  thus,  with  indices  (1,  2,  3)  -B- 

(f?,0,$),25-71 


7r©  —  7r$  —  1/-R)  7ee  —  -R,  7©$  —  cot©, 

p  2  0  (!67) 

7$$  =  —  -Rsin  0,  =  —  sin  0  cos  0. 

Since  the  only  nonvanishing  component  of  D  is  radial  and  since  7pp  =  0,  the 
spray  and  nonlinear  connection  coefficients  derived  from  it  vanish  identically  for 
this  problem: 


GA  =  h bCDBDc  =  \rnnD  ■  D  =  0,  NA  =  GA  =  dBGA  =  0.  (168) 


Denoting  B'  =  dB /dD1  =  dB/dD,  nonzero  components  of  Cartan’s  tensor  in  Eq. 
7  are 


Cm  —  \B' ,  C\ 22  —  C212  —  \B' R" ,  0i33  —  0313  —  \B' R2  sin"  0, 


0221  =  -\B'R\  0331  =  -\B'RA  sin"  0. 


(169) 


The  trace  of  Cartan’s  tensor  in  the  radial  direction  will  be  used  later: 


Cra  =  <7  +  7  +  7  =  3B'/(2B). 


(170) 


The  reference  configuration  space  is  (pseudo)-Finslerian.  Invoking  the  Chern-Rund 
connection  with  vanishing  nonlinear  connection  coefficients  from  Eq.  168,  hori¬ 
zontal  coefficients  are  equal  to  Levi-Civita  coefficients  derived  from  Gab ,  while 
vertical  coefficients  vanish  by  definition  of  the  Chern-Rund  connection: 

BBc  =  =  TfiC  =  7  bc'i  Vbc  =  ^ bc  =  0-  (171) 


The  spatial  configuration  manifold  is  specified  as  (m  :  r  =  x1  G  [0,  r0],  6  =  X 2  G 
[0, 7r],  <f>  =  x3  G  (— 7T,  7r] }.  Metric  tensor  g  necessarily  depends  on  x,  but  by  con¬ 
struction  does  not  depend  on  the  internal  state  vector,  whose  spatial  representation 
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d  is  presumed  radial  but  does  not  further  enter  the  problem  at  hand.  Definitions  and 
identities  of  Section  2.2  result  in 


{r,  9 ,  0}  =  {a;1,  x2 ,  x3},  { d 1,  d2,  d3}  =  {d,  0, 0}.  (172) 


The  following  usual  spherical  form  of  a  Riemannian  metric  tensor  g  LA)  is  assumed: 


g(x) 


1  0  0 

0  r2  0 

0  0  r2  sin2  9 


(173) 


Nonzero  symbols  resulting  from  Eq.  173  are,  with  indices  (1,  2,  3)  (r,  9 , 0), 25,71 

Ire  =  7 '%  =  1  /r,  Yee  =  ~r,  7^  =  cot  9 , 

Yu  =  ~r  sin2  9’  7  u  =  ~  sin  6  cos  9- 

Since  the  metric  is  Riemannian  as  opposed  to  Finslerian  (i.e.,  no  dependence  on 
{d}),  Cartan’s  tensor  vanishes,  nonlinear  connection  coefficients  from  the  spray 
vanish  by  arguments  akin  to  Eq.  168,  and  the  horizontal  Chern-Rund  coefficients 
coincide  with  Eq.  174: 


Cabc  =  0,  <  =  o,  r“c  =  Ybc  •  (175) 

Analogously  to  Eq.  171,  the  following  horizontal  and  vertical  connection  coeffi- 
cents  are  imposed  for  gradients  of  basis  vectors: 

HI  =  Klc  =  r“c  =  Yci  vbc  =  n:  =  0.  (176) 


Motions,  deformations,  and  director  gradients  defined  in  Section  2.3  reduce  as  fol¬ 
lows  under  spherical  expansion/contraction,  with  <p  —  r  denoting  deformation  in 
the  radial  ( R )  direction: 

r  =  ip(R,  D)  =  r(R,  D),  9  =  9,  0  =  <f>;  d  =  d(R,D),  D  =  D(R); 

(177) 


r  dr(R,D)  .  dr(R,D)  dD(R) 
OR  1  dD  dR 

0 

0  ' 

'FYR,D) 

0 

o" 

F(R,  D)  = 

0 

de 

do 

0 

= 

0 

1 

0 

0 

0 

d(j> 
c)<f>  - 

0 

0 

1_ 
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J(R,D)  =  Vg(r,  9)/G(R,  9,  D)F^(R,  D)F2F*  =  [r2FrR(R,  D)]/(R2B3/2)- 

(179) 

C\  i  =  (FR)2,  C'22  =  r2 ,  C*33  =  r2  sin2  0,  Cu  =  C13  =  C23  =  0;  (180) 

det  C  =  J2G]  (181) 

D\x  =  Dyp,  =  d rD  —  NR  +  KRRD  =  dD/dR  =  D' .  (182) 


Internal  state  variable  D  is  physically  identified  with  radial  microscopic  opening  in 
the  material  associated  with  cavitation.  The  defect  associated  with  D  could  be,  for 
example,  a  pore  or  void  in  a  rock  or  metal  or  a  site  vacancy  in  a  crystal  lattice.72,73 
A  lattice  expansion  measure  a,  deformation  gradient  F,  deformation  tensor  C,  and 
Jacobian  determinant  J  are  defined  for  future  use  as  follows: 


a(R,D,D')  =  (183) 


F(R,D,D’)  = 


a(R,  D,  D')  0  0 
0  1  0 
0  0  1 


(184) 


Cab  =  F%gabFB;  Cu  =  a  ,  C22  =  r  ,  C33  =  r  sin  9, 

C12  =  C'i.3  =  C23  =  0; 


(185) 


J  =  yj det  {Cab) /G  =  det  F^J/G  =  ar2/R2.  (186) 

Local  lattice  expansion  occurs  at  R  for  a(R)  >  0,  contraction  for  a(R)  <  0.  Here, 
il}[D(R)]  is  a  continuous  scalar  function  of  its  argument.  As  in  Sections  4  and  5, 
constant  l  with  dimensions  of  length  and  a  normalized  order  parameter  £  are 


£  =  D/l,  r  =  D'jl- 


a  =  5Rcpr  - 


dd  d  D 
dDdR 


Fr  -  wr . 


(187) 


The  constant  l  will  later  be  identified  as  the  value  of  radial  microdisplacement  at 
which  a  material  point  within  the  sphere  supports  no  radial  stress.  Letting  k  denote 
a  constant  depending  on  the  material,  a  more  specific  form  of  the  pseudo-Finsler 
metric  in  Eq.  164  is  henceforth  invoked: 


G{X,D)  =  G(X)B(D)  =  G(X)  exp(2/c£/3),  B(D)  =  exp[(2fcD)/(3/)]; 

(188) 
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B'  /  B  —  diGii/Gu  —  92G22/G22  —  ^3^33/^33  —  2k/ (31).  (189) 

This  can  be  viewed  as  a  Weyl  transformation  or  Weyl  rescaling48  of  the  spherical 
Euclidean  metric  Gab ■  The  radial  component  of  the  trace  of  Cartan’s  tensor  of  Eq. 
170  is 

C%A  =  3B'/(2B)  =  k/l.  (190) 


The  length  of  a  referential  line  element  in  Eq.  14  and  the  corresponding  volume 
form  in  Eq.  15  become 


dX\2  =  exp (2fc£/3) (dR  •  d R  +  R2 d@  •  d@  +  R2  sin2  0d<f>  •  d<f>), 
dQ  =  exjp(k^/3)R2  sin  0d R  A  d@  A  d<f>. 


(191) 


For  £  >  0,  expansion  occurs  when  k  >  0  and  contraction  when  k  <  0.  Physical 
justification  is  obvious  for  the  former:  cavitation  and  void  growth  are  associated 
with  expansion.  The  latter  condition  would  apply  for  collapse  of  existing  pores, 
for  example.  Because  G  is  not  homogeneous  of  degree  zero  in  D,  the  reference 
configuration  space  is  a  pseudo-Finsler  space.16 

6.2  Energy,  Thermodynamic  Forces,  and  Balance  Laws 

Analogoulsy  to  that  considered  in  Section  5.2,  a  compressible  neo-Hookean  strain 
energy  function  W  depending  on  lattice  deformation  tensor  C  is  assumed  to  enter 
^  of  Eq.  52: 


i/(CAB,  Da ,  DA,  Gab )  =  W[CAB(CAB,  DA),  DA ,  GAB ]  +  f[DA,  DA,  Gab], 

(192) 

where  function  /  accounts  for  intrinsic  energy  of  spherical  defects,  such  as  sur¬ 
face  energy  of  a  void  or  cavity.  With  //  and  A  denoting  the  usual  isotropic  elastic 
constants, 


W  =  \{\h(CabGab  —  3)  —  yuln  J+  lA(ln  J)2](l  —  £)2,  (193) 

/  =  (T /l)(DAGABDB/l2  +  Da  GACGBDDfD).  (194) 

A  generally  nonsymmetric  stress  tensor  P  associated  with  W  is 

A*  =  =  \pF^bGAB  +  (A  In  J  -  -  f)2.  (195) 


Approved  for  public  release;  distribution  is  unlimited. 


41 


Since  G(X)  does  not  depend  on  D,  SAB  =  0  follows  from  the  prescription  in  Eq. 
192. 


For  the  present  case  of  spherical  symmetry  with  Eqs.  180  and  185  now  applied, 

d{F^[a(FrR,  D')}}/dFbB  =  (, da/dFrR)SabSA  =  SabSA  =►  PA  =  PA.  (196) 


This  identity,  with  F,1  =  a,  F.j  =  F33  =  1  and  J  defined  in  Eq.  186,  produces  the 
only  nonzero  stress  components 


P  =  Pi  =  PrR 

P'i  =  Pe  =  p'i 


(1-02 

a 


[/r(a2 


pj  =  (1  -  o2 


1)  +  A  In  J] , 


+  A  In  J 


(197) 


The  balance  of  angular  momentum  in  Eq.  53  is  trivially  satisfied  since  F  and  g  are 
of  diagonal  form  and  all  shear  stresses  vanish. 


Using  derivations  from  Section  6.1,  the  following  form  of  Eq.  192  applies  hence¬ 
forth: 


<l,  =  i,(Fl,D,dD/8R)  =  lV[a(^,f),r,e]  +  /«,?') 

(19o) 

=  (i-e)2w„[a(j^,r),r]+/K,e'). 

Adding  /  to  the  elastic  strain  energy  W  reduced  under  spherical  deformation  leads 
to  the  following  reduced  form  of  total  free  energy  density  ip\ 


^(a,  r,  O  =  (1  -  02W"o(a,  r)  +  T^fl  +  W)2 

=  (1  —  ^)2{ ^/z(a2  +  2 r2 /FI2  —  3)  —  n  In (ar2  /  R2)  (199) 

+  lX[ln(ar2/R2)}2}  +  T\e /l  +  lie)2]- 


Nonzero  thermodynamic  forces  are  then  computed  as  follows,  with  stress  compo¬ 
nents  consistent  with  Eq.  197: 


P  =  - — —  [/r(a2  -  1)  +  Aln(ar2/F2)], 

CL 

T  =  Tr(P°  +  P*]  =  -  1)  +  Aln(ar2/fl2)]| 


(200) 


Q  —  Qr 


dip 

~dD 


1  dip 

l~d£ 


(201) 
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z 


d'lp 

Wi 


lchp  1  Of  (l-fl2Wo[a(i%£V]9a 
l  d£'  l  d^’  l  da  d£' 


(202) 

The  following  scalar  function  will  also  be  invoked  later,  where  ^  is  the  function  in 
Eq.  199: 


<a,r,Ul)  =  [2B/(3B')}GABd1GABi;(a,r,Z,e)=2ij(a,r,Z,e) 
=  2(l-02W0(a1r)  +  2T[e/l  +  m2}- 


(203) 


Momentum  balances  specific  to  the  present  problem  are  derived  as  follows.  Thelin- 
ear  momentum  balance  in  Eq.  59  yields  the  three  equations  (a  =  1,  2,  3  =  r,  6,  (f>): 

s.4 P*  +  dBPfdADB  +  PB^B  -  PcA%aFbA  +  P?CcBcdADB  =  0,  (204) 

For  a  =  2,  3,  since  Pf  =  P®  from  Eq.  197,  and  8  =  0,  0  =  <f>  from  spherically 
symmetric  deformation,  these  reduce  to  the  two  trivially  satisfied  equations 

Pftie  -  =  Pf  S  -  Pf  cot  0  =  0,  Pf-faK  =  »■  <205) 


A  nontrivial  linear  momentum  balance  remains  for  the  radial  direction  a  —  1  —  r 
in  Eq.  204: 


dBP  +  WM  +  lp-l(P‘+P^ 


dD 
RA~dR 


+  PCA^  =  0. 


(206) 


The  only  nontrivial  component  of  the  micromomentum  balance  in  Eq.  60  is  for  the 
radial  direction,  C  =  1  =  R: 


az+2z 

m  +  Rz 


/az  aV 

\8D  8D 2 


dD  W 
~dR  +  ~2B 


zm^  =  Q- 


(207) 


Substituting  from  Eqs.  190  and  200  leads  to,  with  P  and  T  radial  and  transverse 
stress  components,  the  following  reduced  form  of  Eq.  206: 


dP 

dP 


R<P-r> 


(208) 
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Substituting  from  Eqs.  190,  201,  202,  and  203  and  multiplying  Eq.  207  by  |  gives 


4  R  4  r  2  \dR  l  de 


=  k 


l-P£+(i-tfw0+je 


(209) 


Relations  in  Eqs.  208  and  209  are  2  coupled  nonlinear  ordinary  differential  equa¬ 
tions  wherein  field  variables  P,  T,  a,  r,  £,  and  Wo  depend  ultimately  on  independent 
variable  R. 


6.3  Problem  Solutions:  Riemannian  Geometry 

Considered  first  is  a  Riemannian  referential  configuration  space  with  k  =  0  — > 
B  =  1  =  constant  =>■  G(X,  D)  — v  G(X).  As  in  Sections  4.3  and  5.3,  kinematics 
are  still  Finslerian  since  <p  can  potentially  depend  on  both  D  and  R  rather  than  just 
R.  Balance  laws  in  Eqs.  208  and  209  reduce  to 


d  P  _ 
dR  ~ 

TIC  + 


n 


2T  /  T 

^'-Te  +  (1-O^o  + 


P 

2~ 


T \  d2r  T 
P )  + 


(210) 

0. 


Unlike  the  problem  considered  in  Section  4.3,  the  first  of  Eq.  210  does  not  neces¬ 
sarily  lead  to  spatially  uniform  stress  fields.  Instead,  ^  =  0  occurs  only  for  the 
case  of  P  =  T,  corresponding  to  a  spherical  and  hydrostatic  stress  state: 


dP/dR  =  0^T  =  P  =  PrR  =  ( R/r)P( f  =  ( R/r)P: f  =  P0  =  constant.  (211) 


From  Eq.  197,  conditions  in  Eq.  211  require  homogeneous  lattice  strain  of  the  form 
a  —  r/R  and/or  complete  cavitation  £(/?)  =  1  \/R  e  S9T  since 

P  =  T  +  (1  -  0V[(«  -  r/R)  -  (1/a  -  R/r )].  (212) 


Solutions  of  Eq.  210  require  more  precise  boundary  conditions.  Two  problems  cor¬ 
responding  to  2  different  sets  of  boundary  conditions  are  addressed:  homogeneous 
damage  of  the  deformed  sphere  over  R  e  [0,  R0\  (i.e.,  microscopic  voids  or  va¬ 
cancies  D  =  Dr  =  D1  evenly  distributed  within  the  volume,  in  conjunction  with 
possible  opening  at  the  origin)  and  localized  cavitation  corresponding  to  a  globally 
stress-free  deformed  state  (i.e.,  complete  cavitation  £  — >•  1  as  R  — >  0). 
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6.3.1  Homogeneous  Damage 

For  homogeneous  damage,  £'(R)  =  0 Vi?  G  [0,  R0]  =$■  £(0)  =  £(Rq)  =  £h-  Bound¬ 
ary  conditions  on  radial  displacement  r  =  p  are  prescribed  as  follows: 

r(0,  D)  =  p0  =  D  =  l£H,  r(Ro,  D)  =  pR  =  aHR0  +  l£H-  (213) 

Here,  pR  is  the  prescribed  outer  radius  of  the  deformed  sphere,  with  aH  =  a(R0 ) 
and  £h  constants.  The  evolving  radius  of  the  discrete  cavity  at  the  origin  can  be 
interpreted  as  po.  The  problem  kinematics,  with  £'  =  0  and  0  =  I),  are  consistent 
with  the  separable  decomposition  of  the  motion  p  into 

p[R,  Z(D)]  =  X{R)  +  #(D)  =  X(R)  +  l£(D),  F  =  FrR  =  dp/dR  =  x!  =  a. 

(214) 

Equations  197  and  210  result  in 

P{R )  =  (1~°^[/r(a2  -  1)  +  Aln(ar2/i?2)]; 

a 

jtH  =  (l-ZH)W0[a(R),r(R)}  +  P 

Given  <pR,  then  it  follows  that  Eq.  212,  the  second  of  Eq.  213,  and  Eq.  215  evaluated 
at  r(R0)  =  <pR  can  be  solved  simultaneously  for  the  homogeneous  damage  field 
£h,  stresses  PR  =  P(Rq)  and  TH  =  T(R0),  and  the  lattice  strain  a(R0 )  =  aR- 
Then,  with  so  determined,  Eqs.  214  and  215  can  be  solved  simultaneously  along 
R  <  Rq  to  determine  the  distribution  of  deformations,  stresses,  and  resulting  energy 
density  inside  the  sphere.  The  total  energy  of  the  elastic  sphere  is  obtained  via  the 
volume  integral 

rRo 

V(<PR,Zh)  =  *h  =  4vr  /  {(1  -  HH)2Wo[a(R),r(R)]  +  T^//}i?2di?.  (216) 

Remark:  Letting  <p(r,D)  — »  ip(R)  =  x(R)  and  modifying  Eq.  213  to  p0  = 
0,  pu  =  ci  hR.o  recovers  a  classical  rather  than  Finslerian  description  of  deformation 
kinematics  in  which  micromotion  D  does  not  contribute  to  macromotion  p.  In  such 
a  case,  the  homogeneous  solution  for  deformation  is  a  —  r/R  —  aR  =  Pr/Rq  = 
constant,  with  P  =  T  =  P0  =  constant  and  £#  determined  by  simultaneous  solu¬ 
tion  of  Eq.  215.  Bifurcation74  to  a  cavitated  state  at  the  origin  is  not  addressed  by 
this  classical  solution. 
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6.3.2  Stress-Free  State 

For  a  stress-free  spherical  manifold  DJt,  P  =  T  =  (NR  e  [0,  Rq],  Particular  bound¬ 
ary  conditions  on  the  order  parameter  £(i?)  or  internal  state  variable  representing 
microdisplacement  are  prescribed  somewhat  analogously  to  Eqs.  98  and  151  as 

£(£oR0)  =  D(e0R0)/l  —  1  —  eQ,  £(-R0)  =  0.  (217) 

As  demonstrated  in  Eq.  218,  the  general  analytical  solution  for  £  is  singular  at 
R  —  0;  therefore,  the  solution  space  is  restricted  to  the  domain  eo  <  R/Rq  < 
1,  where  e0  C  1  is  a  small  parameter.  Linear  momentum  balance  corresponding 
to  the  first  of  Eq.  210  is  trivially  satisfied.  Noting  that  W0  vanishes  for  a  stress- 
free  state,  the  second  of  the  governing  equations  in  210  becomes  the  homogeneous 
nonlinear  second-order  ordinary  differential  equation  for  the  field  £  =  £(f?)  with 
corresponding  general  solution 

£"  +  (2 /R)i'  -  £//2  =  0  =>  £{R)  =  (1/R)[c\  exp(R/l)  +  c2exp(-R/l)\.  (218) 

With  particular  boundary  conditions  of  Eq.  217  imposed,  the  complete  solution  for 
Rq  =  1  is 

««>  -  ei (1  leMm/l)  - eM  (219) 

Remark:  ETnlike  the  problems  considered  in  Sections  4.3.2  and  5.3.2,  the  displace¬ 
ment  condition  at  the  outer  (inner)  boundary  R0  (eoRo)  is  not  arbitrary,  since  van¬ 
ishing  stress  in  Eq.  197  requires  a  —  r/R—  1 WR  :  £  ^  1.  The  boundary  conditions 
in  Eq.  217  thus  imply  r(R0 )  =  R0  in  this  case.  Since  the  radial  order  parameter 
gradient  need  not  vanish,  consistency  with  the  nonlinear  elastic  constitutive  model 
requires  d  =  0  in  the  first  of  Eq.  183  since  a  =  f)n r  =  1  for  this  example. 

6.4  Problem  Solutions:  Finslerian  Geometry 

Considered  now  is  the  general  scenario  wherein  the  referential  configuration  space 
is  pseudo-Finslerian  with  k  ^  0  in  Eq.  188.  Here  the  referential  metric  is  pseudo- 
Finslerian,  and  deformation  kinematics  are  of  Finsler  character  because  motion 
function  r  =  tp(R,  D)  can  potentially  depend  on  D  as  well  as  R.  Balance  laws 
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in  Eqs.  208  and  209  apply  in  full;  substituting  the  first  into  the  second  results  in 


2  TV  T  P 

rie+^'-j(i+km+  2 


21 


R 


T 


P 


+  —f 

de\ 


(220) 


+  (1-O[l-A:(l-O]^o  =  0. 


Solutions  again  require  specification  of  boundary  conditions.  The  same  2  problems 
considered  in  Sections  6.3.1  and  6.3.2  are  now  revisited  in  the  context  of  the  pseudo- 
Finsler  metric  in  Eq.  188,  recalling  that  k  >  0  accounts  for  microscopic  dilatation 
in  the  cavitated  or  porous  zone  omitted  in  the  Riemannian  metrical  representation 
of  Section  6.3. 


6.4.1  Homogeneous  Damage 

For  homogeneous  damage,  £'(R)  =  0 \/R  G  [0,  R0\  =>•  £(0)  =  £(-R0)  =  £h-  Bound¬ 
ary  conditions  on  displacement  are  identical  to  Eq.  213: 


r(0,  D)  =  tp0  =  D  =  l£H,  r(R0 ,  D)  =  (pR  =  aHRo  +  1£h,  (221) 


where  (pR  is  the  prescribed  outer  radius  of  the  deformed  sphere,  with  aH  and 
constants,  and  with  tp0  the  cavity  radius  at  the  core  of  the  body.  Equations  in  214 
still  apply.  Equations  197  and  220  result  in 

P(R)  =  (1  ^(a2  _  1)  +  A  ln(ar2/i?2)] ,  (222) 

CL 

T  PI  (  T\ 

y(l  +  k£H)£H  =  —  f : 1  - -J  +(l-{H)[l-k(l-{H)}Wo[a(R),r(R)l.  (223) 

Given  displacement  boundary  condition  <pR  in  Eq.  221,  the  solution  procedure  is 
analogous  to  that  described  in  Section  6.3.1  (albeit  here  with  k  ^  0),  and  the  total 
energy  of  the  cavitated  elastic  sphere  is  again  obtained  via  the  volume  integral  in 
Eq.  216. 

Shown  in  Fig.  4  are  £  =  £H,  P{Rq)  =  Ph,  and  ^{Rq)  versus  applied  radial  dis¬ 
placement,  computed  via  Eqs.  222  and  223.  Material  parameters  are  identical  to 
those  invoked  in  Section  4.4.1:  //  =  109N/m2,  T  =  lN/m,  l  =  10-9m,  R0  =  103/, 
and  0  <  k  <  In  2.  Trends  are  similar  to  those  of  Fig.  1  and  Fig.  3:  for  fixed  k,  £ 
increases  monotonically  with  increasing  radial  displacement  (Fig.  4a),  radial  stress 
P  increases  to  a  maximum  and  then  decreases  (Fig.  4b),  and  local  energy  density 
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W  at  the  surface  of  the  sphere  increases  monotonically  (Fig.  4c).  As  k  increases,  £ 
tends  to  decrease  for  1.0  <  <Pr/Rq  <  1.5,  while  P  and  w  tend  to  increase.  Varia¬ 
tions  in  fields  with  displacement  are  more  abrupt  in  the  present  spherical  solutions 
than  in  the  1-D  solutions  of  Sections  4.4.1  and  5.4.1.  Peak  radial  stress  and  ap¬ 
plied  displacement  at  which  peak  stress  is  attained  both  increase  significantly  with 
increasing  k,  implying  an  increase  in  cavitation  resistance  and  stability  of  the  ma¬ 
terial  commensurate  with  microscopic  dilatation  represented  by  k  >  0.  Increases 
in  stress  and  energy  correlate  with  a  decrease  in  order  parameter  since  both  P  and 
strain  energy  density  W  are  affected  by  a  multiplication  factor  of  (1  —  £)2  in  Eqs. 
193  and  197. 


Fig.  4  Spherical  deformation,  homogeneous-state  solutions,  l/Ro  =  10  3:  (a)  order  parameter 
£  =  D /l,  (b)  normalized  radial  stress,  and  (c)  normalized  boundary  energy  density 
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Remark:  Letting  (p(r,  D )  — >  </?(i?)  =  x(-R)  and  changing  Eq.  221  to  p0  =  0,pR  = 
cihRo  recovers  a  classical  description  of  kinematics  in  which  micromotion  D  does 
not  contribute  to  macro-motion  (p .  Then  the  homogeneous  solution  for  deformation 
is  a  —  r/R  —  aH  —  Pr/ Ro  =  constant,  with  P  =  T  =  P0  =  constant  and 
determined  by  the  simultaneous  solution  of  Eqs.  222  and  223.  Again,  this  homo¬ 
geneous  solution  does  not  address  the  bifurcation  problem  of  cavitation  as  often 
studied  in  classical  nonlinear  elasticity,74  nor  does  it  address  the  possible  instability 
associated  with  breaking  of  spherical  symmetry.75 

6.4.2  Stress-Free  State 

As  in  Section  6.3.2,  for  a  stress-free  state  P  =  (NR.  e  [0,  R0\.  Boundary  conditions 
on  the  order  parameter  £(R)  are  as  in  Eq.  217: 

£(£oi?o)  =  D(e0R0)/l  =  1  —  eQ,  €(Rq)  =  0.  (224) 

The  linear  momentum  balance  in  Eq.  208  is  trivially  satisfied,  and  ri  =  0  as  re¬ 
marked  in  Section  6.3.2.  The  micromomentum  balance  in  Eq.  220  becomes,  with 
W0  =  0,  the  following  homogeneous  nonlinear  second-order  ordinary  differential 
equation  for  state  variable  field  £  =  £(f?): 

£"  +  (2 /R)?  -  (1//2)(1  +  *;£)£  =  0.  (225) 

This  equation  has  no  known  general  analytical  solution  for  k  ^  0.  Therefore,  solu¬ 
tions  are  obtained  numerically  via  a  second-order  accurate  finite  difference  scheme 
with  an  iterative  evaluation  of  the  nonlinear  term.  Shown  in  Fig.  5  are  profiles 
of  £  computed  via  such  a  scheme,  with  =  0.01,  R.(i  =  10/  and  the  same  range 
of  Weyl  scaling  factor  k  considered  in  Section  4.4.1.  Regardless  of  k,  £  decreases 
rapidly  from  its  maximum  at  R  —  £qRo  with  increasing  R  to  the  imposed  boundary 
value  £(i?0)  =  0.  Increasing  k  provides  no  apparent  change  in  £  for  R  <  R0  for 
spherical  solutions  in  Fig.  5,  in  contrast  to  1-D  Cartesian  solutions  shown  in  Fig.  2. 
The  total  energy  per  unit  surface  area  of  the  sphere  is  shown  in  column  2  of  Table  2, 
where  R.q  =  1  for  normalization.  This  energy  T F  increases  slightly  with  increasing 
k,  similarly  to  trends  observed  for  uniaxial  tension  and  simple  shear  deformations 
in  column  2  of  Table  1 . 
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0.0  0.2  0.4  0.6  0.8  1.0 

R/R0 

(a) 

Fig.  5  Spherical  stress-free  solutions,  l/R0  = 
0.02i?0 


R/Rg 

(b) 

0.1:  (a)  £:  0  <  R  <  R0  and  (b)  £:  0  <  R  < 


Table  2  Stress-free  spherical  solutions  for  l/L0  =  0.1:  total  energy 


Weyl  Scaling  Factor 

Cavitation:  't'i?/(47ri?gT)  •  103 

k  =  0 

1.1077 

k  =  In  | 

1.1080 

k  =  In  | 

1.1082 

k  =  In  l 

1.1083 

k  =  In  2 

1.1085 

7.  Conclusion 

A  new  theory,  in  general  considering  a  deformable  vector  bundle  of  pseudo-Finsler 
character,  has  been  posited,  wherein  the  internal  state  vector  of  pseudo-Finsler 
space  is  associated  with  microdeformation  of  a  material  with  internal  structure.  The 
general  objective  has  been  a  physically  meaningful  theory  that  is  more  descriptive 
and  more  predictive  than  existing  models,  without  ad  hoc  equations  or  numerous 
fitting  parameters.  Rather,  the  focus  has  been  development  of  general,  and  at  times 
more  sophisticated,  governing  equations  instead  of  rudimentary  additions  to  ex¬ 
isting  model  frameworks.  The  specifically  proposed  problem  solutions  offer  new 
physical  insight  into  coupling  of  microscopic  dilatation — captured  herein  by  a  con¬ 
formal  transformation  of  the  metric  tensor — with  fracture  or  slip  in  solids.  It  has 
been  demonstrated  how  the  present  results  can  encompass  known  phase  field  solu- 
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tions  when  a  Riemannian  rather  than  pseudo-Finslerian  metric  is  used  and  how  such 
results  compare  favorably  with  those  of  Griffith’s  fracture  mechanics.  Furthermore, 
the  model  predicts  an  increase  in  peak  strength,  displacement  at  instability,  and 
failure  energy  with  increasing  microscopic  dilatation  in  the  intensely  damaged  or 
localized  shearing  zone,  in  agreement  with  physical  observation.  The  general  the¬ 
ory  is  capable  of  addressing  more  diverse  physical  phenomena  in  condensed  matter 
depending  on  differently  assumed  forms  of  the  fundamental  tensor  (e.g.,  anisotropy 
or  directional  rescaling)  and  different  sets  of  connection  coefficients,  as  may  be 
demonstrated  in  future  work. 

Further  remarks  on  how  the  proposed  Finsler-geometric  framework  may  be  ex¬ 
tended  and  applied  to  problems  of  relevance  in  the  context  of  prior  and  ongo¬ 
ing  work  by  the  author  in  topical  areas  of  defect  mechanics,  structural  transfor¬ 
mations,  and  shock  physics  are  in  order.  Regarding  defect  mechanics,  the  inter¬ 
nal  state  vector  could  be  enlarged  to  account  for  various  components  of  defect 
density  tensors  (e.g.,  dislocations,  disclinations,  and/or  point  defects38'39,72,76-78). 
Phase  field  descriptions  of  deformation  twinning60,79  and  its  competition  with  frac¬ 
ture43,80,81  could  be  modeled  via  identification  of  state  vector  components  with  twin¬ 
ning  shear  as  well  as  local  crack  opening.  Solid-solid  phase  transformations  such 
as  stress-induced  amorphization82,83  are  a  natural  application  of  the  present  frame¬ 
work,  which,  as  has  been  shown,  encompasses  and  extends  existing  phase  field  the¬ 
ory83  to  account  for  additional  physics  such  as  microdilatation.  Particular  crystalline 
materials  of  interest  that  display  twins,  slip/shear  bands,  and  fractures  include  mag¬ 
nesium,60,80  sapphire  or  corundum,80,84  and  boron  carbide.46,82,83  For  shock  physics 
applications,  consideration  of  alternative  nonlinear  elastic  potentials85-87  as  well  as 
inertial  effects  becomes  important.  Explicit  time  dependence  of  field  quantities  has 
been  omitted  in  this  report,  which  has  focused,  for  simplicity/brevity  of  presenta¬ 
tion,  on  an  incremental,  quasi-static  variational  model.  Formulation  of  a  complete 
dynamic  Finsler-geometric  continuum  theory  with  dissipation  poses  no  foreseeable 
difficulties.  For  example,  kinetic  equations  extending  the  Ginzburg-Landau88  or 
Allen-Cahn89  formalism  for  state  vector  evolution  to  the  present  Finsler  modeling 
framework  should  be  readily  possible. 
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